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Introduction 


The  past  five  years  of  this  project  have  been  devoted  to  development 
of  qualitative  theories  of  the  fundamental  electronic  structure  of  pure  and 
impure  semiconductors.  The  basic  theoretical  tools  and  also  conclusions 
relating  to  chiefly  photoemission  type  experiments  and  to  some  extent 
spectroscopy  of  pure  II-VI  and  III-V  compounds  are  presented  in  Chapters  1 
and  2  of  this  report.  Basic  complete  work  relating  to  adsorbates  and  their 
spectra  are  reported  in  Chapter  3.  Several  related  studies  are  being  com¬ 
pleted  at  this  point  and  several  articles  relating  to  these  studies  are 
being  prepared  for  publication. 


\ec\'/ 
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Introduction 


This  investigation  has  the  purpose  of  performing 
self-consistent  energy  band  calculations  on  some  of  the  II-VI 
compounds,  such  as  cadmium  sulfide,  zinc  oxide,  and  zinc  sulfide. 

There  are  many  practical  and  theoretical  reasons 
for  interest  ir.  these  materials. 


and  CdS ) , 
Cu  3  or  C 


These  materials  have  applications  as  phosphors  (ZnS 
as  infrared  detectors  (ZnS),  in  photovoltaic  cells  (cdE- 
iS~Cd7o  he tero June t ior.s  }  ,  ir;  catteries  (ZnO),  as  FETs 


2 


(Cc'S),  in  heteroj 
amplifiers  (CdS), 
They 

well,  as  exam  pier, 
the  covalent  III- 
insulntors.  They 
used  on  the  light 


unction  lasers  (CdS),  and  even  as  acoustic 

arc  of  considerable  theoretical  interest  as 
of  wide  band  gap  .solids,  intermediate  be  Troon 
V  semiconductors  and  the  more  ionic  I-VIJ 
serve  as  a  test  of  calculaticr.nl  methods  previously 
or  II 1-7  compounds  and  ether  compounds  simpler 


V]  cor  cr-u. ids  . 


than  these 


In  particular,  the  occupied  d-levels,  which  lie  in  energy  near  the 
top  of  the  valence  levels,  complicate  the  picture.  The  positions  of 
the  d-bands  have  not  been  correctly  predicted  in  previous  ab  initio 
calculations,  and  in  fact  some  methods  fail  qualitatively  by  not 
obtaining  the  proper  ordering  of  the  valence  s,p  and  d  levels. 

The  II-VI  compounds  have  also  been  of  theoretical  interest 
as  excitonic  systems.  High  densities  of  such  electron-hole  complexes 
have  been  produced  and  observed  in  CdS .  The  electron-hole  liquid  has 
been  observed  and  studied  in  CdS  as  well. 

In  addition,  some  very  interesting  magnet ic  properties 
In  CuCl  and  CdS  have  recently  been  observed.  This  is  the  background 
of  these  experimental  studies.  Anomalous  diamagnetism  has  been 
observed  in  CuCl  under  rapid  pressure  and  temperature  cycling.  In 
some  samples,  up  to  50$  of  the  magnetic  flux  is  excluded.  This  occurs 
at  high  temperatures  ( 1 50  k!  )  and  cannot  be  explained  by  ordinary 
diamagnetism.  The  only  currently  known  phenomenon  that  can  explain 
this  is  superconductivity. 


9 


It's  clear  that  this  would  have  to  be  superconductivity 
of  different  origin.  Acoustic  phonons,  the  known  mechanism, 
are  too  low  in  energy  to  bind  electrons  at  such  high  temperatures. 
Models  using  some  other  form  of  interaction  to  pair  electrons 
seem  to  be  necessary.  Some  form  of  excitonic  binding  is  proposed 
in  many  of  these  models.  Cne  of  the  most  well-known  such  models, 
the  Allender-Bray-F.ardeen  model,  requires  a  metal  and  a 
semiconductor  to  exist  in  very  close  proximity.  Since  CuCl  can 


disproportionate  via  the  reaction  2CuCl 


Cu  +  CuCl  ,  this  model 
2 


has  annealed  to  many.  I-Iowver,  cadmium  sulfide,  which  has  been 
observed  to  have  only  or.e  valence  state  and  is  not  believed  to 
disproportionate,  has  shown  similar  anomalous  diamagnetism.  The 
lav;  of  parsimony  is,  therefore,  against  the  metal-semiconductor 
sandwich  idea.  Also,  in  both  CuCl  and  CdS,  impurities  clearly 
play  a  major  role.  Sufficiently  clean  samples  display  no 
interesting  behaviour  at  all.  And  to  make  life  even  more 
Interesting,  cadmium  sulfide  also  becomes  a  ferromagnet  under  a 


high  (nO  kilogsuss)  applied  field. 


( 

Another  model,  proposed  by  Bishop  and  Overhauser  (  ),  involves 

the  interaction  mediated  by  the  optical  phonons.  This  mechanism 
has  deep  attractive  potential  wells,  about  10  mev  in  depth  for 
a  spacing  of  250  angstroms.  This  suggests  the  possibility  of 
bound  states,  and  the  pairing  is  considerably  stronger,  around 
a  factor  of  ten,  than  the  pairing  from  acoustic  phonons  in 
the  typical  Cooper  pair  in  the  highest  temperature  superconductors 
hno'.-rn,  at  approximately  20  degrees  Kelvin.  This  idea  is  even  more 
intriguing  in  the  light  of  the  fact  that  both  CdS  and  CuCl  are 
strongly  polar  compounds.  This  suggest  that  these  phenomena 
would  not  be  seen  in  the  less  polar  III-V  compounds  and  in 
elemental  semiconductors  such  as  Ge  and  Si.  They’ve  been 
investigated  most  thoroughly  and  indeed,  such  effects  have  not 


beer,  observed. 


Chapter  Cne 


The  calculations  discussed  here  are  based  on 

Hartree-Fock  theory,  the  most  common  approximation  to 

the  exact  non-relativistic  theory  of  a  many-eleotron  system. 

I'll  discuss  the  theory  and  its  applications  to  crystalline 

solids.  The  specialization  of  Hartree-Fock  theory  used 

in  these  calculations,  the  method  of  local  orbitals, 

deserves  and  gets  a  charter  of  its  own. 

The  original  problem  is  finding  exact  solutions 

to  the  Schrodinger  equtation  for  a  nany-electron  system. 

This  is  impossible,  if  the  universe  is  as  small  as  we  think 

it  is.  Consider  the  wavefunction ,  which  is  a  3N  dimensional 

function  if  there  are  N  electrons.  Dividing  each  axis  into, 

say,  100  units,  in  order  to  numerically  integrate  and 

6n 

differentiate,  we  need  to  record  10  entries.  The  same 
problem  arises  if  we  try  to  describe  the  v;avefunction  by 


sets  of  orthogonal  functions. 


6n 

10  is  just  too  big  a  number. 

Quantum  mechanics  is  agaiarst  us.  The  upper  limit  for  the 

2 

rate  of  information  flow  in  a  computer  of  mass  K  is  Me  /h 
(h  s=  Planck's  constant,  c  =  speed  of  light),  in  bits  per 
second.  It's  not  hard  to  see  that  macroscopic  physical 
systems  cannot  b.e  exactly  simulated  by  digital  computers. 

I  emphasize  the  v;ord  digital,  since  the  electron  is  an 
excellent  analog  of  itself. 

Such  a  digits!  simulation  v;ould  be  useless  even 
if  it  could  be  done.  As  suggested,  the  real  solid  simulates 
itself  perfectly.  The  goal  of  physics  is  to  explain  some 
set  of  occurrences  in  terns  of  a  few,  relatively  simple 
ideas,  rather  than  making  a  model  isomorphic  with  the  universe 
and  just  as  confusing. 

We  begin  to  clear  away  some  complications  by 

using  the  Ecrn-Cr.penheimer  approximation.  This  neglects 

any  relationship  between  the  motion  of  the  electrons  and 

the  motion  of  the  nuclei.  This  a  valid  simplification,  since 

3  5 

the  nuclei  range  from  10  to  10  times  heavier  than  an  electron. 


Y/e  can  express  the  /.'-electron  wave  function 


as  a  function  of  the  electronic  coordinates  alone* 


t  =  Y(x,  y  ...  ,  ) 

n  a  1  ^  ^ ;  a 


(1.1) 


Here  contains  the  position  and  spin  coordinates  of  the 


ith  electron.  Y/e  v;rite  the  Schrodinger  equation  for  this 


vavefunction : 

H  s  %(*..*»  ■«„)  (1-2> 
Here,  K  is  the  nor.-relativistic  Hamiltonian  expressed  in 


atonic  units.  In  this  system  of  units,  =  1,  c  ~  2 ,  and 


the  cast  of  the  electron  is  one-half.  The  unit  of  distance 


is  cr.e  bohr  (.529  angstroms)  and  the  unit  of  energy  is  the 


rydberg  (13*6  ev). 

H  -  H  t  H  t/j 

e  n  *nT 


(1.3) 


K  ~  J  l 

h'-l  n  (J^, 

it  '  ' 


(1. 4) 


Ja  X  //?  -  *  / 

*  ~ 

H  independent  of  the  electron  coordinates,  and  has  been 
n 

assumed  be  o  constant,  as  p or  the  Horn-Oppenhc i  mer 


r.  r  :  ro:;i t  ion . 


{ 


The  upper  case  characters  refer  to  nuclear 
properties:  Z  is  the  atomic  number ,  R  is  the  nuclear  position 
(  usually  assumed  to  be  fixed,  like  Pimlico)  and  M  is  the 
nuclear  Eass.  The  lower  case  characters  describe  properties 
of  the  electrons.  y' ,  is  the  coordinate  of  the  ith  elec.ron. 

A. 

The  atonic  numbers  of  the  atoms  under  consideration 
are  relatively  low,  and  since  in  any  event  it  is  the  valence 
electrons  that  are  of  primary  interest,  the  non-relativistlc 
nature  of  the  Hamiltonian  may  be  acceptable. 

The  fundamental  approximation,  r.ov;  applied,  is  the 
independent  particle  model.  In  truth,  in  the  real  solution, 
the  variables  are  not  separable.  We  must  assume  that  they  are, 
or  more  exactly  that  the  true  wave  function  can  be  well- 
approximated  by  this  model.  The  independent  particle  model  says 
that  the  electron  is  acted  on  by  the  average  of  the  other 
electrons.  The  H-electron  wavefunction  is  expressed  as  a  product 
of  cne-electron  wave functions ,  or  as  a  linear  combination  of 
such  products.  A  simple  example  is  the  Hartree  v;avef unction : 


(1.5) 


Since  electrons  are  fermions,  this  [.'-electron  wnvefunctlon  ir.ust 
be  antisymmetric.  Any  interchange  of  two  particles  must  reverse 
the  sign  of  the  wsvef unction.  The  simplest  way  to  antisymmetrize 


is  to  use  the  determinant  of  a  matrix  whose  elements  are 


spinorbitals. 


%  =  (fif'VDc-t 


K)  0/0  ■  ■ '  O/ 
$,(*.)&  I* J  ■:  ' 


(1.6) 


&u  ) 

A/  A‘ 


:;ot  all  such  S^.  are  spin  eigenstates,  but  fortunately  v.*e  have 
here  closed- sho! 1  systems  "hat  car.  bo  expressed  as  single- 
determir.ental  wavof  motions . 

For  any  arbitrary  ,  wc  can  find  the  expectation 

V 

value  of  the  Hamiltonian.  It  is  a  sura  of  one-  and  tvjo-  particle 
3  ntegrals : 

<  V 1,1  V  ~~<1<X  1  fa,l‘ >  h 

I  A  .  -  e*;  6  (>')  <1.8) 

•J  J 

P  is  a  permutation  operator,  constructed  so  that 


r  i  a  >  -  / 


A/  t/ 


,  The  terms  i -- j  are  not  a  problem  since  the 


sol  f-cner/y  cr.r.oo Is  the  • sc  if -exchange •  . 


The  expectation  value  of  the  Hamiltonian,^  V  /  H )  V  \  can  be 

N  As 

minimized  by  varying  the  spinorbitals.  It  has  been  shown  that 

<  j  H  /  %  >  z  E  (i.9) 

<  t  rt> 


where  E  is  the  energy  of  the  exact  ground  state,  .. 

We  vary  the  orbitals,  but  take  care  to  keep  them 
orthogonal,  for  convenience  and  without  loss  of  generality.  This 
is  done  by  using  Lagrange  multipliers.  The  objective  is  to  obtain 
the  lowest  upper  bound  to  the  exact  non-relativistic  energy. 

We  define  the  functional 


L  »  <  t,  >  V  /  %,)  J.  ?yj<  )  d.l 

A  j  S 

and  require  that  the  variation  of  L  =  0  for  1  =  1 , 2 , . .  .  N  . 


0) 


This  yields  a  set  of  H  coupled  nonlinear 


in tegrodif f erent ial  equations  known  as  the  Hartree-F,ock  equations. 

fv*  'J  Spn  +  i  (Jr«j )}  4  M  -/v*ll,1V 


J  and  K  have  the  familiar  Coulomb  and  exchange  operator  forms. 

C/V/  }  0.M 


J. 

**  l  ly-y'i 


a- =/  f  &*(>')#;  o.m 

/  r-y  / 

These  equations  may  be  written  in  matrix  form; 


(1.12) 
(1.13)  ' 


F  0  «  S  £ 


(i.iM 


P  is  the  I.'xN  matrix  of  the  rock  opera  tor .  For  closed 


shells,  it  is  Kermitian,  so  S  may  be  diagonalized  by  a  unitary 
transformation.  Fortunately,  the  11-VI  compounds  have  closed 


shells.  We  may  write,  F  =£./>> 


d.15) 


Combining  this  expression  with  equation  1.7  gives  the  Hartree-Fcck 
ground  state  energy, 

f:'  =  (1/2)  Y  f  €•  .  (1.16) 

KF  W  1  *  J 

The  £>  are  the  single-particle  energies,  almost.  Kooprnan's 

A 

theorem  states  that  given  a  variati orally  stationary  state 

formed  from  N  simultaneous  spinorbitr.l  ( as eucio ) eigenfunctions 

of  the  K  particle  Pock  operator  ?  ,  the  states  Y  and  , 

Aj  f  /  AJ-J 

formed  by  altering  the  list  of  occupied  eigenfunctions  of  F  by 
one  entry,  ore  stationary  with  respect  to  further  variation  of 
that  spinorbital  whose  occupancy  has  been  altered.  If  we 
neglect  relaxa ti on  effects,  £•  is  the  entire  change  in  system 


energy  when  ari  electron  is  added  or  removed. 


In  order  to  equate  the  £~.  $  with  si  ngl  o- part  tele 


energies,  we  must  also  assume  that,  relaxation  effects  are  small  - 


occupied  orbitals  change  only  slightly  when  an 


that  the  other 


There  are  such  changes,  since  the  Pock  operator  is  a  functional 
of  all  of  its  own  solutions,  and  changing  one  of  them  must 
change  all.  The  question  is,  how  much?  These  relaxation  effects 
will  be  discussed  in  a  later  section. 


The  Method  of  Local  Orbitals 

The  method  of  local  orbitals  is  a  variational 

technique  in  which  solutions  to  the  Kartree-Fock  equations 
* 

are  sought  for  a  snail  subsystem  of  the  system  of  interest, 

with  the  caveat  that  all  such  solutions  found  must  be  contained 

within  the  occupied  Fock  space  of  the  original  system.  Sets  of 

such  local  orbitals,  which  are  not  orthogonal,  can  then 

be  assigned  to  each  subunit  of  the  larger  system  using  the 

translational  invariance  of  the  Pock-Dirac  density  matrix, 

thus  spanning  the  occupied  Hilbert  space  of  that  system. 

A  single  noniterative  rotation  of  the  local  orbitals  within 

this  occupied  space  then  yields  the  exact  self-consistent 

solutions  for  the  original  Fock  operator. 

If  £  0f  /(  z  !}  2j  -  'h  is  the  set  of  local 

orbitals  satisfying  the  Hartree-Fock  equations  for  an  n- electron 

system,  then  F  O'  ~  X  /?  (3.1) 

J  * ;  J 

The  Foci: -Dirac  density  operator  for  the  equation  is 


/ 


p  -  5'.  £ 

S  is  the  Inverse  of  the  overlap  matrix  S. 


(3.2) 


j.  -  J  3  />rJ  £<v;  */*  no) 

~  -V 

*  2  *t 

Hho  (y)  is  a  projection  operator  {y  =j>  and^>  =  j>  ) 

chosen  such  thatj>$=  d?'  for  any  orbital  ^'ln  the  occupied  space 

ana yX-f>  for  any  orbital  X  in  t>ne  virtual  space.  This  means 

that  for  any  one-electron  operator  L,  the  projection  of  L 


onto  the  Kartree-Fock  manifold  will  satisfy 


6:  = 


(3.*0 


rie'-p,y^0r  (3-4) 

fo:'  any  occupied  orbital  Oj  ,  and^>  I-y-^-O  for  any  virtual  orbital/. 


Another  approach,  rigorously  developed  by  Gilbert  (  ), 

begins  by  introducing  modified  Hartree-Fock  equations  of  the  form 

(  f  (3.5) 

Above,  A  is  an  as  yet  unspecified  Hermltian  operator  and  rho 
is  the  'diagonal*  (  x  equals  x  prime)  form  of  the  Fock-Dirac 
density  operator  (  ).  We  see  that  rho  is  idempotent,  Kermitian, 

and  projects  onto  occupied  Fock  space.  }3ecause  of  this  last 
property,  it  can  be  shown  that  the  occupied  eigenfunctions  of 


/( 

the  modified  Haa  troo-Foc}:  equation  lie  entirely  v;l  thin  the 
occupied  Fock  space  of  the  unmodified  Kartree-Fock  equations. 

So,  J?  a  ~~  l  0  if  (3.6) 

There  is  no  Kooprr.an's  theorem  for  this  modified 

equation,  and  so  the  £*,  s  cannot  be  interpreted  as  being  nearly 
the  particle  energies,  as  they  can  in  the  regular  Kartree-Fock 
equations.  The  total  energy  of  the  subsystem  still  has  significance 
since  u  unitary  Ft ansf orrn t.1 on  within  occupied  Fock  space 
loaves  the  trace  of  .Fermi tier,  operators  such  as  the  Hamiltonian 

i  V  •  •  -  -.%-*• 

X  -•  •  '  •  v-  . 

The  j  report! os  of  the  // vary  v;ith  the  choice 

A 

of  A.  '.'.'e  tire  loch,  ins  for  functions  that  Eire  localized  on  specific 
lattice  sites  of  crystal,  in  order  to  decide  which  operator 
best  fulfills  such  a  purpose,  we  need  to  consider  the 
min  imization  of  a  functional  [  0  ^  C\t  £?  J 

subject  to  the  variational  constraints  <  //j>  -<.t/  =0; 

p !a)  -!*)  -0;  =o  (3.7) 

J  cf  /' 

can  use  the  method  of  ]/e  mn-.o  multipliers.  The  integral 


constraints  of  equation  3-7  are  taken  into  account  by  the 


parameters  /(  ,  ,  and  the  constraints  that  involve  explicit 
functions  of  X  in  equation  3*7  are  taken  care  of  by  Lagrange 


functions  /£  (x)  and  ?£^{x). 


We  obtain 


({  Ml  ej  •  } 


.2-  (< xfj)  —  < .i )  )  ~2  (j>Dy~)Sy)  (3.8) 

V-  _/g  '*  "/ 

"  /*  X  ""2  2  ,/^>  ^  \  Cj)  <?+>>(*)  ~  O 

K  y- ,  /  /<  ^  S' 


Kultiplyir.g  on  the  left  by  (1-p),  we  see 


0-1 ,2, .  .  ,K 


ar.d  d-jO^*  }  y  *0  for  all 

and  so  j  + /)  (*)  =>  p  (  )  *  $ L  o  J  ~  ?j  (x) 

/f  /£  <y  ^ 

Substituting  equation  3.H  into  3*8,  we  get 


(3.9) 


(3.10) 


(3.11) 


Li 


-  I  2  '/> 


'  *  jy  */- 

•which  may  be  rewritten  as  £ //t>  =  T  ^  /  ■>  > 

t  jb  * 

0-'  J 


wr.ere 


6-  jfcjjA) 


(3.12) 


(3.13) 


(3.14) 


/  unitary  transformation  gives 

0  I  fi’>  =  i'  l  «'> 

Ai 


(3.15) 


Since  we  now  know  how  to  minimize  any  chosen  Hermit lan 


functional  ^  ,  v:e  can  now  got  Lack  to  the  true  problem.  Let 

A  a  G  -  F,  SO  that  £  (  f  hy  A y  &]  ~  O  (3.16) 

We  can  now  minimi z e^lejand  solve  equation  3*5  at  the  same  time, 

and  so  may  pick^Hf-Jto  satisfy  physical  reasoning. 

We  choose  A  so  that  the  solutions  are  localized  on  a  site, 

usually  by  picking  A  to  be  some  potential  well  centered  on  the 

site.  The  Fock  operator  car.  be  separated  into  F  ,  the  atomic 

s 

Fock  operator  for  the  site  s,  and  U  ,  an  external  potential 

s 

operator  for  the  site  s .  we  usually  pick  the  localizing  potential, 

A,  to  be  -U  .  The  local  orbitals  equation  is  now 
s 

f  A  +  V.  -  f  r  3  0  .  "  t  0,  (3.17) 

>  ->  J  J  J  s*  s*  •>> 

This  isn't  an  obvious  improvement  over  the  original 
Hartree-Fock  equation,  but  it  allows  a  systematic  approximation 

that  simplifies  greatly  the  Hartree-Fock  problem  without  inducing 

* 

undue  error.  The  left  side  of  equation  3-17  can  be  analyzed  by 
the  order  of  the  intersite  overlap.  The  Fock~Dirac  density 

(iiy  f/  ~f  y. 

operator  is  p  -■  2  T  0  V  0 
J  s/  7j  Tj, 

with]/  the  iritomtoxic  overlap  matrix. 


(3.18) 


V..  =  f  &*(*)  0.  lx)  r/x 

As  discussed  In  Lowdin  (  ),  V  m*ay  be  expanded  in  a  power 


(3.^9) 


-/ 


Sv 


series  in  V  ;  V  =  3  r)  -  V  +  . . . 

,  -t~  ,  \  •  v 

V  /  <  M 


V>-  5'^ 


(3.20) 


Expending  the  local  orbitals  equation  and  ditcardirig  second 


order  and  greater  in  intersite  overlap  (  )  gives 

IF  +■  a/  -j,  K'pp  1  ^  =  f ' \  Osj 


(3.21) 


and 


,2  1  ^  '> 


where  o  =  )A,  Q  + 

J  ‘r  '  V 

4'«  1  {  '2  ?  _U_ 

5  no  r/,?r-F' 

As  self-consistency  is  approached,  the  last  two  terns  of  equation 
3.21  tend  to  cancel.  Therefore,  the  Fock  space  rotation  changes 


Bloch  wave-functions  into  atonic-like  wavefunctions  that  are 


eigenfunctions  of  the  atomic.  Fork  operfitor. 


We  procede  to  localize  a  physically  realistic 


number  of  electrons  cn  each  site.  In  cadmium  sulfide,  we  choose 

► 

to  have  46  electrons  on  the  cadmium  site  and  18  on  the  sulfur, 
approximating  on  ionic  charge  distribution  (  overlap  and  diffuse 
valence  orbitals  permit  complete  covalency  if  needed,  however). 


We  have  achieved  several  advantages.  The  problem  is 


computationally  simplified,. 


since  the  equation  being  iterated 


does  not  have  to  portray  the  full  complexities  of  the  solid, 
and  since  the  local  orbitals  are  recognizably  modified  atomic 
orbitals  and  aid  in  picturing  chemical  bonding.  This  approach 
has  use  as  a  beginning  in  looRing  at  amorphous  substances,  which 
have  only  a  short-range  order. 

The  local  orbitals  are  calculated  using  the  program 
called  LO?/;S  written  by  hur.z.  A  basis  function  expansion  of 
Slater  orbitals  (SICs)  is  used  for  the  radial  part  of  the  local 


orbital  while  a  spherical  harmonic  describes  the  angular  behavior 

/  (P)  -  Y"(  Od)  I  <r  Ar  .  r*>  1“^K(3.22) 

•<-  3  j  J 

Wc  usually  take  the  /!/  and  jh,  .  from,  vagus  et  al  (  )  where 


optimized  basis  sets  have  been  computed  for  the  atomic  or  ionic 


system  appropriate  to  the  solid  in  question.  In  many  cases,  we 
will  change  this  basis  set  to  achieve  localization,  considering 


at  the  sane  time  the  total  energy  of  the  subsystem. 


The  expansion  coefficient  C  •  is  solved  by  the  matrix 


method  of  booths an  (  ) . 


V.'e  have  only  used  the  long-range  part  of  the  external  potential 

as  the  localizing  potential.  So,  for  ionic  substances,  the 

electrostatic  potential  will  be  screened  while  the  short-range 

effects  of  the  inner  shells  will  change  the  local  orbitals.  Some 

anions  will  not  localize  in  this  situation  since  it  is  the 

fade lung  part  of  the  potential  that  stabilizes  them.  Since  A  is 

essentially  arbitrary,  we  are  free  to  alter  it  in  whatever  way 

i 

that  achieves  optimal  localization,  and  we  often  use- — jl  Jpj 

in  equation  3*21  Instead  of 


CKA  FT  Eli 


o 


Energy  Band  Theory 


The  Hartree-Fock  energy  bands  are  obtained  by  plotting 
the  one-electron  energies  derived  in  equation  1,1.6  against  the 
crystal  momentum  k.  Our  calculations  make  use  of  the  symmetries 
of  the  crystal.  In  the  case  of  these  calculations,  it  is  the 


zinc-blend e  crys to  1  s true ture . 


Bloch's  theorem  requires  that  the  wovefunc tions  of  the 


solid  satisfy 


A  *  *■ 


v?  n  ci'  c 


( y 


R)  = 


AA~  (  >"  ) 
hT 

s*-  (  r) 

h  X 


(2.1) 


for  any  lattice  translation  vector  F;.  The  crystal  momentum  is 

again  k  and  n  is  the  band  index.  V/avefunctions  that  satisfy 

this  condition  are  Bloch  wavefunc tions , 

If  v;c  lock  at  the  set  of  R  ,  we  can  see  that  there 

n 

is  an  element  of  R  that  has  o  minimum  length,  corresponding  to 

n 

the  fundamental  lattice  spacing.  Considering  this,  and  the  fact 
that  R  is  a  group  under  translation,  we  can  sec  that  there  is 


a  rravji;-:"  /  when  v:e  f  ra  n«;  f  or-r  .into  i:  om.  annum  sroce. 


Therefor**,  there  is  b  finite  region  of  non-equivalent  k-polnts  - 
the  first  Erillouin  zone. This  r :«y  be  reduced  to  on  *  irreducible 
wedge '  by  symmetry  considerations. 

The  solution  to  the  h’artree-Fock  equations  is  lengthy 
for  crystals,  and  the  successive  approximation  nature  of  the 
calculations  requires  that  the  equation  be  solved  several  tiroes. 
The  Hartree-Fock  equation  can  be  rewritten  in  terms  of  the. 


Foe V: -Dirac  density  matrix (  ): 


2  |  //f-r/  J  T10J.  r/,'J  &  <x) 

~2J  s* u  )r!/  ~  ^ 

1 t  c{ 

vr,;;rc  J)  ( Sjjt')  ~  ■l  &.(/)  C-\  i >y 


(2.2) 


(2.3) 


If  a  self-consistent  density  matrix  is  available,  then 
equation  2.2  need  only  be  solved  once  to  obtain  the  eigenvalues 
ar.d  eigenfunctions.  Adams  (  )  and  Gilbert(  )  have  shown  that 

the  density  matrix  (x,x’)  is  the  same  for  equations  2.2  and  >.21 


y (/.>-)  -I  =11^  V  £l  (2.*o 

So  we  have  the  required  density  matrix,  and  can  obtain 


crystal  wnvefunctlons  from  the  already  obtained  local  orbitals 
by  constructing  the  rock  operator  and  solving  equation  3.21  just 


or.ee.  The  rat'oix  elements  or  m  are  calculated  to  first  order  in 


inter?, to.uic  overlap  ,  consistent  with  the  calculation  of  the 
local  orbitals. 

The  Bloch  functions  &  (r)  ore  expanded  in  s  basis  set 

KU 

in  which  the  basis  functions  have  the  form 


This  linear  combination  of  basis  functions  technique  is  similar 
tc  the  famous  linear  combination  of  atomic  orbitals  technique 
(LCi-C)  except  for  the  fact  that  the  are  not  free-aton  orbitals, 
but  are  the  local  orbitals  obtained  previously.  Note,  the  local 
orbitals  are  occasionally  enriched  for  the  case  of  virtual  states. 
The  L'FAS  program,  will  not  indicate  basis  functions  appropriate 
for  a  virtual  state  if  that  virtual  state  is  of  a  different 
angular  momentum,  type  than  the  occupied  states;  for  instance, 
e  case  in  which  the  occupied  states  are  all  s  and  p-like,  while 
the  first  virtual  state  is  d-like.  In  this  case,  basis  states'” 
for  the  virtuals  are  added,  consisting  of  spherical  harmonics 
multiplied  by  single  STO's  that  are  chosen  to  have  small  overlap. 

Since  the  crystal  momentum,  k  is  a  good  quantum  number, 
the  lock  sj.ace  in  diagonalized  into  set  a  rate  spaces  for  each 
reciprocal  lattice  vector. 


Vne  internals  necessary  to  perform 


1$ 

the  ca  3  cula  t ion  oro  i:-lnderx?v  hrt  and  need  only  be  performed 
once.  I'.ulticenter  integration  is  done  by  the  Lov.din  c<~f unction 
expansion  method  (  ).  One  of  the  sites  is  chosen  as  center  and 

all  functions  are  expanded  in  terms  of  spherical  harmonics  centered 
on  that  site.  The  calculations  are  performed  using  the  programs 
K-A!3  and  KZCXE  written  by  Kunz.  The  output  consists  of  energies 
ar.d  coefficients  of  the  basis  functions  at  20  selected  points 
of  the  Brlllouin  zone. 


.26 

•  Correlation  Corrections 

As  previously  mentioned,  Rartree-Fock  theory  does 
not  give  exact  solutions  to  the  true  many-electron  Hamiltonian. 

Since  we  are  using  a  single  determinental  wavefunction ,  electrons 
are  affected  only  by  the  mean  field.  Electrons  of  the  same  spin 
have  some  of  their  true  pair  interactions  taken  into  account, 
but  there  is  no  pair  correlation  at  all  btween  two  electrons 
of  opposite  spin.  It  is  obvious  that  the  Coulomb  force  between 
two  electrons  is  independent-  of  spin,  and  should  keep  electrons 
from  the  near  neighborhood  o*'  any  other  electron. 

V,'e  say  that  the  notion  of  electrons  of  opposite 
spins  is  uncorrelated,  and  that  electrons  ’with  like  spins  are 
incompletely  correlated.  The  correlation  energy  Is  usually, 
defined  as  the  difference  between  the  energy  obtained  from 
our  rr.esn-field  Rartrco-Fock  theory  and  the  exact  non-relatlvistic 
energy  of  the  system.  This  exact  energy  could  in  principle  be 
determined  by  using  a  computational  technique  called  configuration 
interaction 


In  this  method,  a  wave function  is  used  that  is  a  linear 
combination  of  determinants.  It's  impractical  for  all  but  the 
smallest  molecular  and  atomic  systems,  and  is  far  too  time- 
consuming  to  be  useful  in  solids. 

For  solids  in  general,  and  specifically  in  the  case 
of  the  II-VI  compounds  under  discussion,  we  attempt  to  estimate 
the  correlation  correction  to  the  Hartree-Fock  calculation.  The 
correlation  correction  for  the  one-electron  eigenvalue  is  the 
error  from  Koopmoru  *  theorem.  (  )  and  from  the  independent 

particle  model. 

According  to  Koopmans *  theorem,  an  orbital  eigenvalue 
is  approximately  the  difference  beyween  the  Kartree-Fock  energy 
of  the  K-electron  system  with  the  level  occupied  and  the  Kartree- 
Fcck  energy  of  the  N-l  electron  system  with  that  same  level 


unoccupied.  These  Kartree-Fock  excitation  energies  are 


-  t  -  E 

ii  nr  t/r- 


(4.1) 


for  states  that  are  occupied  in  the  Kartree-Fock  ground  state 


t.  rlA/} 

r  -  [  -  F 

\-c'  jf  *  h  r- 


(4.2) 


for  virtual  stater  of  the  Karfree-Foc.k  ground  state. 


When  Wfi  attempt  to  improve  the  Kartrec-Fock 


ZT 


approximation,  we  v.ish  to  keep  the  general  band  scheme,  so 


(*')  _  (/!(-•) 

€  ^  t  ~  t 

Y\k. 

rW) 


(4.3) 


(4.4) 


where  the  eigenvalues  are  the  exact,  eigenvalues  of  the  system. 
vie  assume  that  the  correlation  energies  ore  a  small  perturbation 


of  the  original  Uartree-Fock  bends,  and  we  can  then  write  E  as 


(L)  (L)  (I,) 

E  =  B  +  E 


(4.5) 


here,  E  is  the  total  correlation  energy  of  a  system  with  L 
c 

electrons .  Equations  4.3  and  4.4  nay  now  be  written  as 

o  <::)  <::-}) 

e  rs  e  +  (K  -  E  )  (4.6) 

nk  nk  c  c 

0  (N)  (k-1) 

e  -  e  +  (E  -  E  )  (4.?) 

nk  nk  c  c 

A  useful  approximation,  valid  for  nonmetals  with 

valence  band  width  Irss  than  the  Hartree-Fock  optical  band  gap, 

«» 

was  developed  by  i'antelides  et  al  (  ).  They  showed  that  equations 


4.6  ar.d  4.7  nay  be  replaced  b$ 

0  (N-1  ) 

K  «  E  +  E  (h) 
nk  nk  nk 


(4.8) 


(4.9) 


In  these  equations,  ,  E  t  (e)  is  the  self-energy  of  on  electron 

nk 

that  occupies  the  one-electron  Kartree-Fock  state  ri*k  in  an 

(N-l ) 

K-electron  system.  E  (h)  is  the  total  energy  change  in  the 

nk 

remaining  N-l  electrons  when  the  electron  occupying  the  state  in 

(N-l) 

question  is  removed.  E  (h)  is  the  self-energy  of  a  hole. 

nk 

Pantelides  et  al  {  )  have  produced  some  model- 

independent  results  concerning  these  self-energies.  They  find 
that  the  self-energy  of  holes  is  always  positive,  and  therefore 
the  valence  hands  from  the  original  Hartree-Fock  calculation 
always  move  up  on  the  energy  scale  upon  correlation.  These  self- 
er.ergies  grow  as  v;e  approach  the  bottom  of  the  valence  bands ,  so 
on  ta lance  the  valence  bands  are  narrowed.  The  self-energies  cf 
the  lower  electrons  in  the  conduction  band  are  negative,  so  the 
conduction  bands  move  down.  The  models  used  have  only  small 
change  in  the  amount  of  shift  for  the  low-lying  conduction  band, 
and  so  we  use  a  rigid  downward  shift  for  these.  We  end  up 
shifting  the  valence  bands  higher  and  rigidly  dropping  the 
conduction  bands.  This  produces  a  smaller  optical  gap  than  that 
oh  twined  by  the  Hartrec-Fock  calculations. 


We  more  nearly  approach  the  experimental  gap  by  this  n9.rrov;lng. 

These  self-energies  were  first  calculated  by 
Toyozawa  (  ).  The  theory  of  this  calculation,  the  electronlc- 

polaron  method,  war.  considerably  further  developed  by  A,  B.  Kunz, 
so  that  it  now  predicts  hole  self-energies  as  well. 

In  the  original  Hartree-Fock  theory  particles  respond 
only  to  the  average  position  of  the  other  electrons  and  ions. 

This  is  obviously  incorrect,  or  at  least  Incomplete ;  the 
independent  charge  v;.i  11  polarize  its  surroundings,  especially  if 


those  surroundings  have  sufficient 


f ,i ms  to  respond.  The  electronic- 


polar  on  model,  dreoso^^  the  coxio  uo  uj on  bar  d  electrons  ana  valence 


band  holes  with  quanta 


of  the  polarization  field.  These  quanta  are 


excitons.  In  the  model  we  are  using,  the  excited  states  of  the 

(N) 

crystal  are  simulated  by  a  di spersionless  band  of  excitons.  E  (e) 

Uk 

(N-l  ) 

and  E  (h)  are  the  interaction  energies  of  a  bare  electron  and 
nk 

hole  v.'ith  this  field.  These  energies  are  called  polarization 


energies.  This  model  uses  second-order  perturbation  theory  to 


calculate  the  self-energi no . 


/ 


V 


►  . 


These  polarization  energies  ore  calculated  on  the 
basis  of  a  model  in  which  the  hole  or  electron  is  fixed  in  space. 
This  method,  the  Kott-Llttleton  method,  uses  a  perturbative 
approach  to  find  the  induced  dipole  moments  on  all  ions  of  the 
crystal  due  to  the  localized  charges.  This  calculation  takes  into 
account  the  field  from  .the  induced  dipoles  -  it  is  self-consistent. 
Since  it  takes  a  finite  time  for  the  crystal  to  respond  to  such  a 
change,  moving  charges  should  induce  less  polarization  than  the 
model  static  charge.  Therefore,  this  calculation  should  give  an 
upper  bound  for  the  polarization  energies  of  actual,  mobile 
electrons  and  holes. 

This  electronic-polaron  model  is  basically  a  long- 
range  scheme.  Only  in  the  limit  of  large  distance  can  v;o  find  these 
changes  by  assuming  that  they  can  be  described  by  dipoles. 

For  short  distances,  the  shape  of  the  wavefunctions  plays  a  major 
role,  and  quadrupole  and  higher  multipole  effects  cannot  be 

* 

ignored.  The  changes  in  the  central  atom  and  its  near  neighbors 
certainly  cannot  be  modeled  by  a  dipole  very  successful ly! 


Short-range  correlation  calculates  the  change  of  the 

nearby  orbitals  when  an  electron  is  added,  or  removed.  To  find  these 

corrections,  v:e  do  simple  atomic  calculations  of  different 

ionization  states,  supplemented  by  cluster  calculations  in  the  case 

of  negative  ions  v/here  the  added  electron  extends  over  a  significant 

region  of  space.  In  the  limit  of  zero  overlap  betv/een  atoms,  the 

short  range  correlr. tion  correction  to  the  hole  energy  is 
sr 

E  (h)  =  -C-Z  (4.10) 

s  c  f 


where  <f  is  the  energy  of  the  level  the  electron  was  removed  from 

and  K  is  the  Rartree-Fock  ionization  energy.  The  ionization 
scf 

energy  is  the  difference  between  the  Hartree-Fock  energy  of  the 

system  with  the  level  occxipiod  and  the  Hartree-Fock  energy  with 

that  same  level  full.  The  short-range  energies  so  calculated 

are  again  upper  bounds,  since  electrons  and  holes  are  not  as 

localized  as  this  model  po trays  them. 

The  Kartree-Fock  calculation  and  the  following 

correlation  corrections  give  a  set  of  one-electron  energy  bands. 
* 

They  must  now  be  compared  to  experiment  to  test  the  adequacy  of 


our  approximation. 


CADMIUM  SULFID2 


c 


At  room  temperature  and  standard  pressure  CtiS 
tnllizes  in  the  zinc-blende  lattice  with  a  lattice  constant 


of  5 • 31 8  angstroms  (10.99  atomic  units)  and  in  the  hexagonal 
fcrr.  with  lattice  constants  of  4,1343  and  6.7234  angstroms  ( 
7 . 31 36  and  12.705^  atomic  units).  The  zinc-blende  form,  the 


subject  of  these  calculations,  is  composed  of  two  interpenetrating 

++ 

face-centered  cubic  sublattices  occupied  by  Cd  and  S  ions, 
displaced  relative  to  each  other  by  1/4  of  the  diagonal  of  the 


unit  cube.  The  symmetry  properties  of  the  zinc-blende  lattice 
2 

(  st-jce  grout)  T  )  have  been  discussed  by  Pnrmcnter  (  ),  and  in 

d 


following  discusr.li 


r ouc-'acrt ,  Smoluchowr.ki  and  Wigncr  (  )  will  be  used. 


We  first  examine  the  energy  states  at  the  prillouin 


10 

zone  center,  the  point.  Cadmium  ions  contribute  a  filled  (4d) 

shell  to  the  valence  energy  region,  and  sulfur  ions  contribute 
2  6 

a  (3s)  (3p)  configuration.  The  p  and  d  levels  hybridize  to 

some  extent,  since  they  are  separated  by  less  than  nine  ev;  still, 

the  valence  bands  are  not  strongly  hybridized.  The  top  valence 

band  is  predominantly  p-like  (  ever  75/0  nnd  the  next  lower  band 

is  predominantly  d-type,  again  over  75/0  We  expect  the  lowest 

conduction  band  to  be  derived  from  the  cadmium  5s  levels,  as  is 

usual  in  compounds  that  have  any  ionic  character.  In  the  zinc- 

blende  structure,  the  crystal  field  splits  the  fivefold  degenerate 

d  levels  into  a  trirly  degenerate  state  and  a  doubly  degenerate 

1 «; 

state.  The  threefold  degenerate  p  states  stay  degenerate  under 

12 

this  crystal  field  and  transform  like  .  The  zinc-blende  lattice 

15 

has  no  Inversion  symmetry,  so  the  bands  at  the  garnma  point  need 
not  have  a  definite  parity.  The  conduction  band  contains 

•5 

both  p  and  d  contributions,  in  fact. 

Several  previous  calculations  have  been  cone.  The  valence 
c  a  1  ] y  the  neper  ones,  ami  the  lower  conduction  bands 


bar 


esr" 


have  been  we? 1-desor ibed  by  various  pseudopotentio 1  calculations. 
That  is,  they  accord  with  cur:  ent  exp^ "er,  since 
the  pseudopotential  method  is  basically  a  parametrization  scheme, 
in  which  the  parameters  are  determined  by  experiment,  little 
critical  new  information  can  be  obtained.  If  the  experimental 
evidence  has  been  misunderstood  or  misinterpreted,  the  •  v 
pseudopotential  will  simply  predict  the  mistakes  or  misunder¬ 
standings  that  it  sprang  from.  A  pseudopotential  calculation 
is  an  aid  to  undei standing,  but  it  tends  to  not  be  falsifiable. 
Such  calculations,  such  as  those  by  Cohen  and  Eergstresser  (  ) 

also  depend  or.  the  validity  of  the  cancellation  theorem,  which 
is  not  exact  and  which  is  much  less  valuable  for  systems  that 
have  localized  states  that  ore  closely  comparable  in  energy  with 
the  valence  states  -  here,  the  cadmium  d  states  are  the  ones 
not  easily  described  because  of  their  local  nature. 

A  related  but  more  theoretically  rigorous  technique, 
the  orthogonalized  plane  wave  method,  has  some  of  the  same  errors. 
In  this  method,  the  valence  and  conduction  states  ore  described 


by  plane  waves  orthogonalized  to  the  core  states,  taken  as  constant. 


in  these  previous  OF,-,  calculations,  as  performed  by  t£uv<ewn  and 
Stukel  (  ),  Euwerna  and  Collins  «t  al  (  )  and  by  Stukel  et  al  (  ) 

there  have  been  errors  of  up  to  three  ev  in  the  p-like  valence 
and  conduction  bands  in  compounds  made  of  first-row  atoms. 

Also,  the  4d  states  are  misplaced  by  about  a  rydberg  -  probably 
related  to  their  core-like  nature.  The  level  ordering  of  the 
cadmium  bd  and  the  sulfur  s  bands  seems  to  be  reversed  from 
that  measured  by  photoemission  by  about  6  ev,  Altogether  it 
is  easy  to  see  that  there  is  a  need  for  a  first-principles 
calculation  that  correctly  predicts  the  major  features  of 
the  band  structures  of  these  materials. 

Therefore,  the  first  all-electron,  self-consistent  '' 
noner.pirical  band  calculation  using  nonlocal  exchange  has  been 
performed  on  CdS.  The  emthods  of  calculation  are  as  discussed 
earlier  ir.  this  work.  It  should  be  mentioned  that  a  selfr-consistent 


band  calculation  has  been  performed  by  Zunger  and  Freeman  (  ) 

using  a  local  density  approximation  for  the  exchange  potential. 


Their  results,  as  well  as  earlier  ones,  will  be  discussed  in  the 


( 

TAHi  .il  u . : 

I 


STC 

BASIS  FUKCTICI' 

!S  FOE  cad: 

a uk  a: 

!D  LOCAL 

OR  DUAL  CO 

EFFICIENTS 

s  Easls 

p 

Ba  s  i  s 

d 

Basis 

n 

z 

n 

Z 

n 

Z 

Oj 

Oj 

1  j 

1  j 

2  j 

2  j 

1 

1 

70.66 

3 

31  .29 

3 

20.99 

2 

1 

45.69 

2 

1  9.62 

3 

11.48 

3 

2 

38.00 

3 

10.45 

3 

7.67 

4 

2 

20.71 

3 

6.83 

4 

5.75 

5 

3 

IS.  45 

4 

7.00 

4 

2.50 

6 

3 

10.05 

4 

4.55 

5 

2. 00 

7 

4 

5.31 

5 

2.90 

3 

5 

2.90 

5 

2.00 

Q 

S 

5 

1.90 

10 

20 

30 

40 

50 

21 

j 

r» 

c 

C 

C 

C 

C 

J 

j 

j 

j 

J 

1 

0 . 04^36 

- . 01 447 

0.05145 

0.03371 

0.0033? 

-.17167 

2 

0!  951 73 

0.37037 

-.C0305 

-.01190 

-.00465 

-.86 02? 

j 

-.03876 

0 . 1 7?49 

0.4:2  31 

0.19171 

0.03936 

0.10571 

1, 

0 . 01 953 

-1 .17556 

-1 . 2i  1  96 

-.54694 

- . l 0?6? 

-.23043 

r 

J 

-.00903 

-.02754 

0.7293? 

0.35372 

0.07401 

0.17254 

s 

r\ 

0.00142 

- . 00827 

0.93'-:50 

0.53579 
-1  .21  004 

0.03914 

-.02615 

7 

-.00065 

0.00355 

-.16742 

-.22539 

0.00925 

O 

c.  00033 

-.00172 

O.05351 

-.0346? 

0.04950 

-.00394 

9  v 

-.00015 

0.00093 

- .  0451  9 

-.00954 

0.95921 

31 

41 

51 

32 

42 

j 

c 

c 

c" 

C 

r> 

J 

j 

j 

J 

J 

- 

i 

-.05253 

o.c:  8i  2 

0.00473 

0.08,700 

- . 03?49 

O 

-.46795 

0 . 1 5690 

0. 01238 

0. 66066 

-.16153 

3 

1 .21641 

-.52352 

C. 05904 

0. 3G001 

-.23294 

u 

-.30972 

-.16943 

-.37952 

0.01590 

0.6524  c 

5 

0.30198 

0./l3795 

0.44531 

-.00499 

0.73322 

( 

-.051  43 

0.78304 

-.12?46 

0.0c 302 

-.15556 

7 

0.0:794 

0.08963 

0 . 3  35  4  3 

1- 

-.00775 

- . 01 56b 

- .  0 

TABLE 


SIC  BASIS  FUNCTIONS  FOR  SULFUR  AX  5  LOCAL  Ci.BI 1AL  COEFFICIENTS 


s  Basis 


p  Basis 


J 

n 

z 

n 

z 

Oj 

Oj 

1  j 

Ij 

1 

1 

17.60 

2 

1  8.00 

2 

2 

15.45 

2 

13.42 

3 

3 

•  9.68 

2 

10.0c 

I. 

2 

5.36 

2 

7.75 

< 

3 

2.59 

2 

4.80 

6 

3 

1.63 

3 

2.32 

7 

3 

1.32 

1  0 

2C 

30 

40  21 

31 

V* 

c 

n 

c 

C  C 

C 

i 

j 

j 

j  j 

J 

* 

0.^228 

-.23268  0. 

062? 0  C 

.02395  0.0004: 

- ■ 04392 

* 

c.:  r-"5: 

- !  1  i s  c  9  0! 

0397?  0 

(>24  ’.i  C  .  02  9  v  " 

0.23470 

*4 

> 

-. coc2a 

0.13650  -. 

02455  0 

.10551  0.001  Vj 

-•  •  53346 

4 

0. DC 2 35 

0.93736 

33199  - 

,30626  0. 33621 

0.57490 

-.CG054 

0.00602  c. 

51 001  0 

."-1224  O.6692: 

-.06509 

- 

C. CCC26 

-.00035  0. 

59220  - 

.62170  0.C0595 

-.I93OO 

? 

41 

32 

-.00074 

-. R4226 

A 

r> 

c 

WJ 

J 

: 

C. 04330 

0.00 

2 

-.22559 

0.00 

3 

0.54462 

0.00 

t : 

-.50114 

0.00 

r 

J 

0.273c? 

■  .  00 

6 

-.4503? 

7 

0.3nC?6 

rolnrizfibll  ltiee 


:d  \y  Tessnan  et  al  (  )  were  then  used  for 


Cd  and  S  along  with  the  dielectric.  :  -^.nt  of  6.32 

to  calculate  the  polarization  energies  using  i.:-  i-.<»tt-Littleton 

method.  We  obtained  values  of  0.1866  ry  and  0.2584  ry  were 

ob-t-a-ired  as  the  polarization  energies  associated  with  the  S 
++ 

and  Cd  holes  respectively.  Being  less  tightly  bound,  the 

sulfur  anions  polarize  to  a  greater  extent  around  a  hole  at  the 
~+ 

Cd  cation  site  than  do  cadmium  ions  around  a  sulfur  hole, 

producing  a  larger  polarization  energy  for  the  cation.  We  also 

calculated  short-range  relaxation  energies  for  the  states  of 

sr  — 

interest  in  the  valence  region;  these  are  E  (S  )  =  0.07824  ry, 

3p 

sr  sr  ++ 

E  (S  )  =  0.C63C2  ry ,  and  E  (Cd  )  -  0.116  ry.  These  corrections 
3C;  4s 

were  a.dded  to  the  Hamiltonian  matrix  which  was  then  rediagonalized 
to  give  the  correlated  valence  bands.  Conduction  bands  were 
computed  by  shifting  the  Hartree-Fock  conduction  bands  by  -.2584  ry, 
the  polarization  energy  of  a  conduction  band  electron  on  a 

+  -T 

Cd  site. 

Shown  in  figure  4.1  are  the  correlated,  energy  bands 
of  Cut. .  The  calculated  band  streture  shows  cadmium  sulfide  to  be 


Figure  4.1  Correlated  energy  bands  of  CdS  for  the  normal 


lattice  constant  of  5*181  angstroms  (  10.994  au). 


a  direct  b'nd  pep  semiconductor  with  the  gao  fit  the  fr» r*:r?isi  point. 

This  conclusion  is  in  agreement  with  the  previous  band  calculations. 
Three  major  non-overlapping  regions  constitute  the  valence  bands 
in  this  system.  The  S  3s  derived  band  lies  19  e'v  below  the 

valence  band  edge  and  is  about  1 .4  eV  wide.  The  next  region 

+  + 

is  primarily  derived  from  the  Cd  4d ,  and  lies  12  eV  below  the 

top  of  the  valence  band.  Tills  band  is  about  one  eV  wide.  The 

uppermost  valence  band  is  primarily  derived  from  5  3P  levels, 

ar.d  is  about  3*9  eV  wide.  The  lowest  conduction  band  is  s-like 

++ 


ar.d  is  derived  from  the  Cd.  5s. 

S 

3s  ami  S 

4  s 

levels . 

The  band  gap  is  found 

to 

be  direct 

and 

equal  to 

7.1  eV. 

The  optical  value  of  the  gap  ,  E 

2.53  eV  ( 

). 

cbta  ir.ed 

from 

experiment,  is  in  serious  disagreement.  The  correlation  model 
here  used,. the  electronic-rolaron  model  and  it  limit, 
fcott-I.  i  ttleton  theory,  is  adapted  to  insulators  and  ignores 
short-range  polarization  effects.  Such  effects  should  be  small 
lr.  atoms  and  in  systems  where  the  local  orbitals  are  only  slightly 
perturbed  froo  the  free  atom  or  i  or. ,  but  in  polar  semiconductors 


f’.: oh  •:  s  C;i-  ,  thor 


guarantee  that  the  local  orbitals  are 


is  no 


Figure  4.3  XPS  spectrum  for  CdS  according  to  Ley  et  al  (  ). 

I  and  I  describe  fine  structure  of  the  highest  valence-band 
1  2 

peak.  S  is  a  shoulder  on  the  high-binding  side  of  this  peak. 

1 


II  represents  the  second  valence-band  peak. 


( 

that  similar  to  the  free  ions.  In  this  class  of  compounds,  there 
is  substantial  charge  density  in  the  internucl ear  region. 

It  has  been  argued  that  such  effects  act  to  narrow  the  band  gap. 

The  density  of  states  for  CdS  as  calculated  from  the 
original  band  structure  calculation  is  given  in  figure  4.2.  For 
comparison,  figure  4.3  shows  the  corrected  XPS  spectrum  obtained 
by  Ley  et  al  (  ).  The  experimental  evidence  clearly  shows  the 

two  peaks  in  the  density  of  states  of  the  upper  valence  band. 

The  positions  of  these  peaks  are  also  in  reasonable  agreement 
with  experiment .  The  upper  peak  lies  1.6  cV  below  the  top  of 
the  valence  band  by  Ley's  measurements,  while  our  calculation 
gives  a  peak  at  approximately  1.4  eV.  Experimentally,  the  second 
peak  is  at  4.1  eV,  while  this  calculation  has  a  peak  at  eV. 

It  also  seems  possible  to  identify  the  shoulder  of  the  upper  peak. 

•» 

From  our  band  calculation,  this  shoulder  seems  to  be  at  about 
2.C  eV,  while  Ley"s  measurements  put  it  at  2.1  eV. 

The  position  of  the  d-like  levels  is  also  correctly 
predicted  3r.  these  calculations.  These  levels,  primarily  formed 
from  cadmium  4d  states,  are  found  to  peak  at9;&4  eV  belov; 
the  top  of  the  valence  band  experlmentrillly ,  according  to  Ley. 


Figure  4.5  Imaginary  part  of  the  dielectric  constant  of  CdS 

(  this  calculation) 


Figure  4.6  Experimentally  derived  imaginary  part  of  dielectric 
constant  in  CdS,  from  Cardona  et  si  (  ). 

Derived  from  ultraviolet  reflection  measurements. 


Cur  calculation  gives  a  value  of  10.F  c7  below  the  tor  of  the 
valence  band,  in  very  good  agreement  with  experiment.  This  is 
in  sharp  contrast  with  previous  pseudopotential  and  OF'./ 
calculations,  which  misplaced  this  band  as  much  as  12  eV  below 
its  actual  position.  It  should  be  mentioned  that  Ley's 
photoeraission  studies  were  done  on  CdS  in  the  hexagonal  form, 
but  since  the  first  two  shells  around  the  central  atom  ore 
Identical  in  the  hexagonal  and  zinc-blende  structure,  the  atom 
is  exposed  to  a  very  similar  potential,  and  no  substantial  change 
in  the  density  of  stn tea  is  expected. 

Ir:  figure  4.4  the  joint  density  of  states  as  calculated 
from  our  band  structure  is  shown.  In  figure  4.5,  the  imaginary 
part  of  the  dielectric  function  is  calculated  from  the  same 
theoretical  band  structure.  In  figure  4.6,  Cardona  et  al  extract 
the  imaginary  port  of  the  dielectric  constant  from  ultraviolet 
reflection  measurements  on  zinc-blende  CdS  (  ).  There  is  clear 

qualitative  agreement,  especially  jlth  the  central  peak  and  the 
two  R’.i n— tv-*i v s  approximately  3  and  5  eV  above  that  central  peak. 

The  toil  be  lew  the  central  pr  «k  lr.  energy  is  exaggerated  due  to 
tie’  eve  res'.  *.  :•  a  te  of  the  band  gap  in  our  calculations. 


t 

To  sum  up:  these  calculations  correctly  predict  the 
cajor  features  of  the  energy  bands  in  CdS,  as  determined  from 
optical  and  other  measurements.  Insofar,  this  first  all-electron, 
self-consistent  non-empirical  method  has  succeeded.  The  d-llke 
levels  are  correctly  predicted,  as  v?as  not  the  case  with  earlier 
calculations.  Nothing  in  the  pure  bulk  energy  bands  suggests 
8n  explanation  of  the  anomalous  diamagnetism,  and  it  would  seem 
that  the  effect  is  not  an  intrinsic  one.  It  may  be  hoped  that 
this  better  understanding  of  the  electronic  structure  of  the 
cure  solid  will  aid  in  the  understanding  of  the  anomalous 
dia~a.rr.etic  state,  perhaps  by  acting  as  necessary  first  step  in 
the  examination  of  the  properties  of  defects  and  impurities  in 
cadmium  sulfide. 


ZINC  OXIDE 


At  room  temperature  and  standard  pressure  ZnO 
crystallizes  in  the  hexagonal  (  wurtzite)  form  with  lattice  constan 
of  3.2^9  angstroms  (6.430  atomic,  units)  and  $.19%  angstroms  (9.S24 
atomic  units).  The  zinc-blende  form  with  identical  nearest-neighbor 
distance  has  a  lattice  constant,  of  6,595  angstroms  (3.634  atomic 
units)  and  is  the  subject  of  these  calculations.  It  is  composed  of 
two  interpenetrating  face-centered  cubic  sublattices  occupied  by 

-i-J_  —  — 

Zn  and  0  ior.s,  displaced  relative  to  each  other  by  1/4  of  the 
diagonal  of  the  unit  cube. 

V.'e  begin  again  by  examining  the  energy  states  at  the 

10 

Erillouln  zone  center.  Zinc  ior.s  contribute  a  filled  (3<i) 

shell  to  the  valence  energy  region,  and  oxygen  ions  contribute  a 
2  6 

(2s)  (2p)  corf iguration.  In  this  crystal  structure,  the  threefold 
degenerate  r  states  stay  degenerate  and  transform  like 


15 


The  crystal  field  splits  the  fivefold  degenerate  dlevels  into  a 

triplv  degenerate  '  state  and  a  doubly  degorern  te  }'  state. 

15  12 

The  zinc  3d  levels  and  the  oxygen  2p  levels  lie  quite  close, 
and  this  calculation  shows  shows  all  eight  valence  bands  lying 
within  a  five  eV  region.  The  zinc-blende  lattice  has  no  inversion 
symmetry,  so  the  bands  at  the  gamma  point  need  not  have  any 
definite  parity.  This  calculation  shovis  that  this  upper  valence 
region  exhibits  strong  p-d  hybridization.  The  band  lying 
beneath  this  upper  valence  region  lies  nearly  24  eV  below  the 
top  c ^  the  valence  band,  and  is  primarily  derived  from  oxygen 
2s  svstes.. 

In  general,  in  ionic  compounds  such  as  ZnO  we  would  expect 
the  lowest  conduction  band  to  be  primarily  formed  from  the 
zinc  4s  level.  The  actual  calculation  shows  that  although  the 
zinc  4s  states  play  a  major  role,  oxygen  2s  ana  js  states  actually 
play  a  larger  role  in  this  first  conduction  band.  The  next 
cond\:ction .  band,  triply  degenerate  at  the  ".one  center,  are  formed 
almost  entirely  from  zinc  4p  states  .This  is  suggested  merely 
by  co" raring  zinc  with  its  successor  in  the  periodic  table,  gallium. 


Several  previous  calculations  have  been  performed  on 
zinc  oxide.  Attempts  using  local  pseudopotential  theory,  such  as 
those  by  Rossler  (  )  and  Bloom  and  Ortenburger  (  )  have  not 

beer,  entirely  successful.  Since  the  pseudoootentinl  theory  depends 
upon  a  cancellation  cf  the  strong  core  part  of  the  potential  by 
the  usual  requirement  that  the  valence  electrons  be  orthogonal 
to  the  core  electrons,  first-rov?  elements  would  seem  to  be 
unsuitable  for  this  approach.  The  point  is  that  valence  p-states 
for  first  row  elements  are  not  required  to  be  orthogonal  to  any 
p-core  states.  Investigators  have  attempted  to  alleviate  this 
problem  by  empirical  adjustments  to  the  pseudopotentiol s ,  but 
it  hasn't  worked  well.  When  nonlocal  pseudopotentials  were  used, 
as  in  the  calculations  of  Pnelikowsky  (  ),  a  better  agreement 

with  experiment  is  achieved.  Even  so,  much  of  the  value  of  the 
pseudcrotential  approach  is  lost.  There  is  not  the  same  confidence 
that  the  pseudopotential  will  retain  its  pred  j  tct.ive  powers  in 
different  compounds,  because  the  physical  rationale  is  weakened. 
The  mad  5  fled  and  nor.  local  calculations  still  suffer  from  all  the 


original  weakness  of  pseudopotert ials :  reliance  on  experiment. 
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Therefore,  we  perform  on  all-electron,  self-consistent 
noner.pirical  band  colculatlon  on  ZnO.  The  methods  of  calculation 
ore  the  same  as'  used  previously  on  CdS. 

The  calculation  was  begun  by  calculating  local  orbitals 
for  7.  inc  oxide  with  the  normal  lattice  constant.  The  zinc  and 
oxygen  basis  sets  of  Bagus  et  al  (  )  were  used  and  were  modified 

for  this  calculation.  Two  additional  diffuse  STOs  were  added  to  the 

s  basis  on  zinc.  The  results  of  the  local  orbitals  calculation 
are  listed  in  tables  5*1  and  5«2  for  zinc  and  oxygen  ions, 
respectively. 

Kartree-Fock  bands  were  calculated  for  20  k  points 
in  the  irreducible  wedge  of  the  Brillouln  zone.  The  usual 
exaggerated  band  gap  appears;  we  proceed  to  apply  the  correlation 
methods  previously  discussed. 

Folarizabill ties  listed  by  Tesstran  et  al  (  )  were*  then 

4* -4-  —  ~~ 

used  for  Zn  end  0  along  with  the  optical  dielectric  constant 
cf  •’.'.036  to  calculate  the  polarization  energies  using  the  Kott- 

Li ttletori  method.  We  obtained  values  of  0.2172  ry  and  0.2864  ry 

— —  ++ 

ss  the  polarization  energies  associated  with  the  0  and  Zn 
holes  respectively .  'We  also  calculated  short-range  relaxation 


igure  5-1  Correlated  energy  bands  of  ZnO  for  the  normal 


lattice  constant  of  4.595  angstroms  (  8.634  au 


energies  for  the  state  of  interest  in  the  v.-> fence  region: 


thiS 


sr  ++ 

ds  E  (Zn  )  =•  0.3!402  ry.  These  corrections  were  added  to  the 

3d 

Hamiltonian  matrix  which  was  then  rediagonalized  to  give  the 


correlated  valence  bands.  Conduction  bands  were  computed  by  shifting 


the  Hartree'Fock  conduction  bands  by  -.2864  ry,  the  polarization 

++ 

energy  of  a  conduction  band  electron  on  a  Zn  site. 


Shown  in  figure  5*1  are  the  correlated  energy  bands 


of  Zr.O.  The  calculated  band  structure  shows  zinc  oxide  to  be 


a  direct  band  gap  semiconductor  with  the  gap  at  the  gamma  point. 


This  conclusion  is  in  agreement  with  previous  band  calculations. 

Two  major  ncn-overlapplng  regions  constitute  the  valence  bands  in 

this  system.  The  highest  region  is  about  4  eV  in  width,  and  is 

divided  into  two  subregions ' with  almost  no  overlap.  The  higher 

end  wider  of  these  subregions  is  ,  approximately  3  eV  in  width, 

— —  ++ 

is  composed  of  O  2p  and  Zn  3d  levels,  while  the  lower  subfcand 

++ 

is  almost  dispersionless  and  is  of  almost  pure  Zn  3d  character. 

The  lowest  valence  band  is  about  2.3  eV  in  width  and  is  composed 

— ~  ++ 

of  0  2p  and  Zn  3d  levels. 


C  ^  r  o  #, 


The  band 
the  gap  is  ex 


gao  is  found  to  be  direct  and  equal  to  1C. 86  eV 
perim.ciitally  found  to  be  3.3  eV ,  it  must  be 


< 


I 
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Figure  $.2  Density  of  states  of  ZnO 
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Figure  5.3  Corrected  XPS  spectrum  for  Zr,0  according 


to  Ley  et  el( 


that  the  conduction  band  structure  here  predicted  is  seriously 
in  error.  Previously  mentioned  potential  problems  in  our  correlation 
model  may  explain  this  discrepancy,  but  it  seems  most  likely  that 
there  is  some  flau  in  the  estimate  of  the  polarization  energies. 

The  density  of  states  for  ZnO  as  calculated  from  our 
band  structure  is  given  in  figure  5*2.  For  comparison,  figure  5*3 
shows  the  corrected  experimental  density  of  states,  derived 
from  x-ray  photoemission  experiments  conducted  by  Ley  et  al  (  ). 

The  d-like  levels  have  been  subtracted  out  from  the  experimental 
data,  since  they  dominate  the  spectra.  The  double  peak  in  the 
density  of  states  of  the  uppermost  band  is  clearly  shown  in 
experiment  and  in  our  calculation.  The  d-like  levels  in  our 
/  calculation  have  some  structure,  are  not  smoothed  into  a  single 

.peak  as  in  the  experiment.  The  experimental  peak  is  centered 

■ 

8. Si  eV  below  the  top  of  the  valence  band,  according  to  Ley, 

Our  d-complex  Is  situated  approximately  5  eV  below  the  top, 
but  the  qualitative  picture,  that  of  a  double-peaked  valence 
band  with  a  very  sharp  d-band  about  3  eV  below ,  corresponds 


closely  with  experiment. 


Our  calculations  also  (rive  a  band  derived  from  0  2s 
states  (not  shown  on  figure  5*2)  centered  around  24  eV  below 
the  top  of  the  valence  band  -  this  corresponds  to  the  lowest 
peak  in  figure  5*3 

In  addition,  we  have  calculated  the  joint  density  of 
states  (  figure  and  the  imaginary  part  of  the  dielectric 


constant  (  figure  5*5)  for  zinc  oxide. 


_£ure  5.^  Joint  Density  of  States  of  ZnO 
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Figure  5-5"  Imaginary  part  of  the  dielectric  constant  of  ZnO 


CONCLUSIONS 


We've  used  ab  initio  Hartree-Fock  theory,  as  well  as 

relaxation  and  polarization  correlation  corrections  to  that  theory, 

to  calculate  the  electronic  band  structures  of  cadmium  sulfide  and 

zinc  oxide.  In  the  calculation  of  cadmium  sulfide,  hybridization 
++  — 

of  the  Cd  ~s ,  S  js  and  S  Us  is  a  major  factor  in  the  lowest 
conduction  band.  The  calculation  shows  that  CdS  is  a  direct  gap 
semiconductor  with  a  gap  of  7.1  eV  at  the  gamma  point.  This  tends 
to  show  that  the  mechanism  proposed  by  Abrikosov  to  explain  the 
enomalous  diamagnetism  originally  seen  in  CuCl  (  and  discredited 
there  by  Weidman  (  )  )  is  not  applicable  to  CdS,  since  it  requires 

a  small  indirect  gap. 

The  calculated  one-electron  energy  bands  are  compared 
with  the  published  optical  data.  The  valence  bands  are  in  excellent 
agreement  with  photoernission  and  reflection  data,  and,  in  particular 


the 


position  of  the  d-like  band  is  correctly  predicted,  in  contrast 
with  previous  calculations. 

The  calculation  for  zinc  oxide  shows  that  there  is 
very  significant  nixing  in  the  upper  valence  bands  between  the 
C  2p  and  Zn  3d  levels.  It  seems  clear  that  previous  pseudopotential 
calculations  that  could  not  correctly  take  into  account  the  zinc 
d  levels  were  incapable  of  explains  the  valence  bands.  The  calculation 
gives  a  band  gap  of  1.0.86  ev,  and  predicts  ZnO  to  be  a  direct  gap 
serai  conductor  with  the  gap  at  the  gamma  point.  In  the  case  of  both 
CdS  and  ZnC ,  it  should  be  recalled  that  we  are  essentially  just  solving 
Dyson’s  equation,  and  that  higher  order  diagrams  are  quite  capable 
of  accounting  for  the  difference  between  the  (exaggerated)  gaps 
from  our  calculations  and  the  measured  experimental  gaps.  The  fact 
remains  that  these  calculations  predict  the  valence  structure? 
the  one  electron  bands  for  zinc  oxide  also  accord  well  with  published 
optical  experiments,  /..gain,  the  position  of  the  d-like  bands  is 
correctly  predicted. 

The  anomalous  diamagnetism  observed  in  CdS  remains 
unexplained.  It  scans  to  be  an  extrinsic  phenomenon  . 
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The  observation  of  ferromagnetism  upon  an  applied 
field,  as  well  as  the  anomalous  diamagnetism,  is  intriguing. 

It  suggests  that  the  pairing  mechnaism  may  favor  a  state  with  a 
spin  of  one,  as  opposed  to  the  Cooper  pair,  with  a  net  spin  of  zero. 
Each  pair  would  have  an  intrinsic  magnetic  moment.  This  possibility 
has  been  discussed  for  classical  phonon-mediated  superconductivity 
by  P.  Anderson  and  F.  Korel(  )  but  seems  not  to  have  been 
observed,  except  possibly  at  extremely  low  temperatures  (  ). 

An  exclton-medlated  pairing  might  well  favor  the  1=1 

or  hisrhar  state.  The  superfluid  state  in  He  is  suggestive. 

3 

As  yet  the  details  of  the  interactions  are  not  know, 
although  it  does  seem  that  chemical  impurities  play  a  role.  It  may 
still  be  possible  to  test  some  of  these  ideas  by  a  phenomenological 
theory  like  the  Glnzburg-Landau  theory.  One  possibility  seems 
interesting;  in  the  Ginzburg  -Landau  theory,  two  characteristic 
lengths  appear,  the  coherence  length  and  the  penetration  length. 

In  such  a  model  revised  for  e  p-wnve  system,  a  third  length  arises, 
the  characteristic  length  for  the  change  in  spin  direction. 

The  two  lengths  of  Ginzburg-Landau  theory  allow  the 
existence  of  a  distinctive  surface,  which  may  have  a  lower  energy 


r/\ 
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than  the  bulk  state,  which  leads  to  the  Interesting  and 
technologically  useful  type  II  superconductors.  This  still  exists 
with  three  characteristic  ler.rths,  but  another  effect  becomes  possible. 

A  second  surface  layer  is  introduced;  we  may  compare  this  to  the 
earth,  with  crust,  mantle  and  core.  If  the  mantle  is  energetically 
favored,  bubbles  of  a  certain  size  would  be  energetically  favored. 

These  suDerconducting  domains  might  explain  the  very  high  but  finite 
conductivity  seen  in  the  anomalous  states  of  CuCl  and  CdS,  and 
alignment  of  these  domains  might  explain  the  feiromagnetlsm  seen 
in  CdS  under  high  applied  field.  At  the  moment,  this  is  all  speculation. 


We  await  further  experiments. 
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ABSTRACT 

Energy  bands  and  ionization  energies  for  electrons  in  A1P,  GaP,  and 
GaAs  are  obtained  using  a  new  method  proposed  by  Kuna,  et  al.'  The  valence 
bands  and  ionization  energies  obtained  are  found  to  be  in  good  agreement 
with  experiment.  The  conduction  bands  vary  substantially  from  experiment. 
Finally,  a  method  for  obtaining  better  conduction  bands  is  proposed. 


1 .  INTRODUCTION 


A  new  ab  initio  method  for  obtaining  correlated  energy  bands  was 
recently  proposed  by  Kur.z,  et  al . 1  The  need  for  a  new  method  becomes 
apparent  v:hen  the  inherent  problems  associated  with  current  ab  initio 
methods  are  considered.  Existing  methods  can  be  broken  into  two  classes, 
those  which  use  the  nonlocal  Hartree-Fock  (H-F)  exchange  potential  and 
those  which  use  some  local  density  approximation  to  the  exchange  or 
exchange-correlaticn  potential.  The  former  class  of  methods  requires 
calculations  which  are  both  lengthy  and  complicated.  Once  such  a  H-F 
calculation  has  beer,  performed,  there  is  still  a  need  to  obtain  correlation 
corrections.  These  corrections  are  poorly  understood  for  covalent 
semiconductors  like  the  III-V  compounds.  On  the  other  hand,  local  density 
approximation  calculations  give  results  which  are  highly  dependent  on  the 
choice  of  the  exchange-correlation  potential.  No  single  potential  has  been 
found  to  give  good  results  for  all  classes  of  compos.  ,s.  For  example, 


Herman,  e_t  a_l_.  2 ,  found  that  in  using  both  the  Slater^  and  the  Kohn-Snam.1' 
exchange  potentials  for  several  semiconductors ,  the  Slater  potential  agreed 
more  closely  with  experiment  for  II-VI  semiconductors,  while  Kohn-Shan 
exchange  gave  better  results  for  the  III-V  compounds.  Ever,  in  cases  where 
a  given  potential  yields  the  correct  valence  and  conduction  band 
structures,  such  calculations  have  not  always  been  able  to  place  the  core 
levels  in  their  correct  locations  relative  to  the  valence  bands.  Finally, 
local  density  calculations  have  not  been  successful  in  obtaining  the 
ionization  energies  for  electrons  in  most  compounds.  We  believe  that  many 
of  these  difficulties  are  a  result  of  using  a  Hartrce  potential  which 
includes  a  sel f-rcpulsion  term,  on  the  assumption  that  the  potential  from 
on»-  electron  in  an  infinite  crystal  will  be  negligible.  As  has  been  shewn 
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by  Kunz,  et.  al.  \  this  assumption  is  not  correct  in  many  cases  and  is  a 
source  of  signi ficar.t  error  for  insulators. 

The  new  method,  which  will  be  referred  to  as  the  Hartree-plus  method, 
uses  the  correct  Hartree  potential,  i..e.  with  the  self-repulsion  removed, 
along  with  a  local  exchange-correlation  potential.  The  Hartree-plus  method 
was  applied  by  Kunz,  et  al . 1 ,  to  the  solid  rare  gases  and  NaCl  yielding 
bands  which  are  in  good  agreement  with  both  experiment  and  previous  H-F 
plus  correlation  calculations.  The  III-V  semiconductors  provide  a  good 
test  of  the  range  of  applicability  for  the  Hartree-plus  method  since  they 
are  wide  band  semiconductors,  as  opposed  to  the  rare  gases  and  NaCl  which 
are  narrovr  band  insulators. 

In  Sec. 2,  H-F  theory  is  developed  and  correlation  is  discussed.  The 
Hartree-plus  theory  and  details  of  the  calculations  are  presented  in  Sec. 3 
and  SecJl  respectively .  Ir.  Sec. 5,  the  results  of  Hartree-plus  band 
calculations  for  A1P,  CaP ,  and  Ca/.s  are  presented  and  compared  to 
experiment  and  previous  theoretical  calculations. 

2.  BASIC  THF, C:-':- 


The  system  of  interest  contains  n  electrons  and  N  nuclei.  Using  the 
Born-Opponheimer  approximation \  and  neglecting  relativistic  effects,  the 


Hamiltonian  is: 


i”  / 


*  i  Z 

I  v/  -  » 


(!) 


where: 


/;  -  -  V;  -2  = 


#/ 


The  energy  is  in  Hyd bergs ,  e^= 2,  uppercase  letters  refer  to  nuclei, 

lowercase  letters  denote  electrons,  and  the  prime  on  the  sun  indicates  that 


the  self-repulsion  term,  i=j,  ir.  excluded.  The  problem  now  ir.  to  solve  the 


Pape  0  2. 


Schrodinger  equation: 

**  sr'Wix.'-x*')  (V 

where  is  the  exact  n-elcctron  wavefur.ction  for  the  i^11  excited  state, 
and  the  are  taken  to  include  both  space  and  spin  coordinates. 

To  simplify  this  problem,  we  approximate  the  exact  wavefur.ction  "ifc <tj 
with  a  single  Slater  determinant  formed  from  a  set  of  or.e-electrcr; 
orbitals,  (^(xj)),  i=1,*”,n  and  j=1,***,n: 

a  ■'■*>.)  •=  (»/)"* Jet  1} &(*,)}! 

The  expectation  value  of  for  the  trial  function  "T*  ,  will  be  a 

rigorous  upper  bound  to  the  exact  ground  state  energy  if  the  resulting 


energy  Eft'1'  is  stable  against  variation  in  the  Performing  such  a 

variation,  subject  to  the  constraint  that  the  be  or thonornal ,  yields  an 
equation  for  the  ft;  : 

rft:(?,)  ~  f  A.y 
where  F  is  the-  Fcck  operator: 


r  -  f,  -r /A,  f v.  ^  Pte,  Ol 

p  (  a ,  xT.  >  -  ^  o?,  >  ir/j  ’ 

P(pt)ft;(/,)  ■=-  ft;  (k3) 


13) 


In  cq.(3),  P(2,l)  is  an  exchange  operator,  the  nay  be  used  as 

operators,  and  is  the  Fock-Dirac  density  matrix.  The  second  tern 

is  kr.ov.-n  as  the  direct  or  ilartree  term,  and  the  last  term  is  the  exchange 
term.  f.’ote  that  the  prime  is  no  longer  necessary  since  the  Kartree 
solf-er.crgy  is  exactly  cancelled  by  the  sel f-exchange.  Diagonalizing  the 
matrix  Aj-j  gives  the  standard  Hartrec-Fock  equation: 

F  ftp  ft)  -  e?P;  17,)  CV; 

The  eigenvalues  of  ec .  ( '• )  are  given  meaning  by  Koopnans 1  theorem?. 
Labelling  occupied  states  with  and  virtual  states  with  2=n+1,***t 


Pace  gr  J 


the  eigenvalues  are: 

c  /  —  tv  tl/ 

__  r=".^t-l)0 _  pi''.'8 


where  E*/  is  the  energy  of  the  single  Slater  determinant  approximation  to 

the  total  m-electron  wavefur.ction  obtained  by  adding  or  removing  the  oc^h 

orbital  to  the  n  orbitals  used  to  form  the  H-F  ground  state.  Thus,  and 
6 

6.  are  the  energies  that  an  electron  would  have  if  its  instantaneous 
location  were  independent  of  the  instantaneous  locations  of  the  other 
electrons,  and  if  no  relaxation  were  allowed.  To  improve  on  the  H-F 
eigenvalues  while  retaining  the  independent  orbitals,  one  defines 
correlated  energies  by  replacing  the  approximate  energies  in  eq.(5)  with 
the  exact  energies: 

(V 


~  _ l*rf) 

wc*  ’*"  *-„•?  tZo 

Tr.e  energies  €;  and  €c,  are  quasiparticle  excitation  energies.  The 
quasiparticles  associated  with  these  energies  are  referred  to  as  holes  and 
conduction  electrons  respectively. 

There  are  two  basic  approaches  to  finding  the  quasiparticle  energies. 
The  better  understood  of  the  two  is  the  H-F  plus  correlation  method,  which 
is  essentially  a  perturbative  approach.  In  this  method,  eq.(4)  is  solved 
first,  and  then  the  correlation  effects  are  added  in  as  a  small  number  of 
corrections.  This  approach  to  obtaining  correlated  bands  was  reviewed  and 
developed  in  a  formal  way  by  Pantelides,  et  al.7  The  correlation 
corrections  are  fairly  well  understood  for  both  metals  and  narrow-band 
wide-gap  materials.  Unfortunately  the  III-V  compounds  do  not  fall  into 
either  of  these  classes.  The  other  approach  is  to  use  an  equation, 
analogous  to  the  H-F  equation,  which  somehow  incorporates  correlation  into 
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the  potential.  This  is  the  approach  taken  by  the  Hartree-plus  method. 

3.  THE  HARTHEE-PLUS  METHOD 

In  the  past,  band  calculations  using  local  exchange-correlation 
potentials  have  solved  an  equation  of  the  form: 

Os.r.  $;(*')  =  (K)  (?) 


where: 

Oif.  ~  V,  +  ( Vac  J 

is  a  single  particle  operator  which  replaces  the  nonlocal  exchange  in  the 
Fock  operator  with  a  function,  Vxc  ,  of  the  local  density ,/CKr) .  Unlike 
the  Fock  operator,  eq.(7)  no  longer  has  an  exact  cancellation  of  the 
self-repulsion  by  a  self-exchange.  It  has  generally  been  assumed  that  the 
self-repulsion  is  negligible  since  it  comes  from  one  electron  in  an 
infinite  crystal.  This  is  in  fact  true  for  an  electron  in  a  totally 
delocalized  Plcch  state.  However,  what  is  actually  desired  is  the  energy 
required  to  add  one  quasiparticle  to  the  ground  state  system.  For  all 
insulators,  and  many  semiconductors,  these  quasiparticle  orbitals  are 
local.  In  this  case  it  is  reasonable  to  work  in  the  local  Heitlcr-Lcndcr. 
representation.  For  such  local  orbitals,  the  self-energy  is  not  in  fact 
small  and  should  be  removed  explicitly. 

In  the  local  representation,  let  $«(r)  be  the  i^h  local  orbital  at 
site  M.  Tne  corrected  potential  then  is: 


where: 


H*-  V,  Vh  (?/  )  +  VKc  if&Q 

Vj  (?, )  -  jVxv  <y,>  pc *3  xl )  ~ 


is  the  nonlocal  Hartree  potential.  Although  the  Hartree  potential  is 
nonlocal,  the  nonlocality  is  easy  to  deal  with  by  using  a  single  Hartree 
potential  for  all  orbitals  of  the  form: 

Y„(W )  -  *./» CX,K>  -r,  /4„>  z?(s.£)<4~I  (?) 
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where: 

Zy*  (s.a)  =•  2  i  ]<?,*■ C'3.)l  ■*  /  )}x  I i£  ■ -  l -jr?A  ef?. 

Here  the  bra-ket  notation  is  used  to  emphasize  the  operator  nature  of  the 
self-repulsion  term. 

For  Vxc  ue  have  used  the  local  exchange  potential  of  Kohn-Shan1*  and 
the  dielectric  screening  function  of  Robinson,  Bassani,  Knox,  and 
SchrieffertriBXS)8.  Thus: 

y«.  f m sc«;  a  o) 

where: 

vrcv  =  -,t[(-£>)/o<?)Ta 

PM = /-  -'(£■) Vw 

and: 

oc  646  EpCF)~]~>A: 

Although  it  might  seem  inconsistent  not  to  remove  the  self-exchange  in 
eq.(ll),  one  can  argue  that  the  self-exchange  will  only  be  a  significant 
fraction  of  the  total  exchange  potential  in  regions  of  low  density  where 
F(&c)  approaches  zero. 

1J.  THE  CALCULATIONS 

To  calculate  the  Hartree-plus  bands  shown  in  Figs.  1-3,  we  first 
performed  an  LCAO  calculation  using  H-F  atomic  wavefunctions  obtained  by 
Bagus,  Gilbert,  and  Roothaan^,  In  this  calculation  Nesbet's  symmetry  and 
equivalence  restrictions^  were  applied.  Also,  the  crystal  potential  was 
approximated  by  the  first  term  of  its  expansion  in  Yim»  i-e.  1=0.  The 
next  tern  in  the  potential  would  be  1=3  and  should  not  have  a  large  effect 
on  the  s,  p,  and  d  orbitals,  which  are  all  that  are  of  interest  here.  The 
potential  and  matrix  elements  are  evaluated  using  10  shells  about  the 
central  site.  V.'e  represent  the  remainder  of  the  crystal  by  a  residual 
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Kadclung  potential.  To  obtain  this  potential,  we  used  effective  charges 
given  by  Phillips11.  This  LCAO  calculation  can  not  be  expected  to  give 
good  results  for  the  conduction  levels  due  to  the  use  of  Hesbet's 
restrictions1®.  The  self-energies  of  the  atomic  orbitals  have  been  removed 
for  both  the  occupied  and  the  unoccupied  parts  of  the  orbitals.  This  means 
that  the  conduction  energies  have  been  obtained  as  though  the  conduction 
electrons  felt  the  potential  from  n-1  electrons,  instead  of  all  n  valence 
electrons.  To  remedy  this  the  self-energies  are  explicitly  reinserted  for 
the  virtual  states.  At  this  point  the  Hartree-plus  matrix  is  no  longer 
diagonal.  Hov;ever  the  coupling  between  the  occupied  and  virtual  states  is 
very  weak  and  a  second  diagonal ization  produces  new  valence  bands  which, 
for  GaP,  differ  in  energy  from  the  previous  bands  by  at  most  Taking 
all  this  into  consideration,  we  believe  that  the  final  Hartree-plus 
operator  so  obtained  is  consistent  with  the  configuration  used.  Finally  we 
form  ?.  Hartree-plus  matrix  using  atonic  orbitals  for  the  core  states  and  89 
planewaves  for  the  conduction  and  valence  states,  and  diagonalize  it  to 
obtain  the  fine]  band  structure.  For  A1P,  we  repeated  this  calculation 
using  27,  51*  and  65  planewaves. 

5.  RESULTS 

In  all  the  band  calculations  performed  here,  energy  levels  were 
obtained  at  21  nonequivalent  points  in  the  Brillouin  zone  for  an  f.c.c. 
lattice.  The  bands  are  drawn  along  axes  of  high  symmetry  which  connect  the 
points  of  high  symmetry  in  the  Brillouin  zone.  The  symmetry  labels  used 
are  consistent  with  Parnenter1^.  In  all  cases  the  origin  has  been  chosen 
to  be  the  cation  site.  In  addition,  the  orbitals  and  their  energies  are 
superscripted  to  indicate  whether  they  are  valence  (v)  or  conduction  (c) 


states. 
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The  valence  band  str  uctures  for  AlP^,  obtained  using  27,  51,  65,  and  89 
planewaves^  are  summarized  in  Table  1.  Wo  experimental  or  theoreti.cal 
values  were  available  for  comparison.  The  bands  obtained  with  £9 
planewaves  are  shown  in  Fig.1.  By  comparing  the  results  for  several 
numbers  of  planewaves,  we  can  estimate  the  size  of  the  errors  due  to 
incomplete  convergence.  For  the  top  two  p-derived  bands,  our  results 
appear  to  be  converged  to  within  0.1  eV.  The  third  p-like  band  is 
converged  to  within  about  0.3  eV,  and  the  separation  of  the  s-band  and  the 
p-bands  is  converged  to  within  roughly  OJi  eV.  On  the  otherhand,  the  basic 
structure  is  stable,  as  is  the  width  of  the  s-band.  We  assume  that  this 
convergence  information  is  also  valid  for  our  other  calculations. 

In  Tables  2a  and  2b,  the  valence  band  structures  for  GaP  and  GaAs 
shown  in  Figs. 2  and  3  are  compared  to  experiments  and  previous  theoretical 
calculations.  It  is  particularly  interesting  to  compare  the  Hartree-plus 
results  with  the  tv:o  first  prir.ciples  0?W  calculations.  For  both  GaP  and 
GaAs,  the  Hartrec-plus  method  gives  valence  band  widths  which  are  in  better 
agreement  with  experiment  than  those  obtained  in  the  first  principles  OPV.' 
calculations,  ar.d  are  comparable  in  quality  to  the  empirically  adjusted 
calculations,  especially  when  the  lack  of  convergence  is  taken  into 
consideration.  All  of  the  other  calculations  give  better  results  than  this 
work  for  the  separation  of  the  s-band  from  the  p-bands.  This  error  in  the 
Hartree-plus  results,  1.0  to  2.0  eV,  is  probably  due  to  incomplete 
convergence  and  relaxation  effects  that  have  not  been  included.  None  of 
the  other  calculations  listed  give  the  position  of  the  d-bands  relative  to 
the  top  of  the  valence  bands.  For  the  Hartree-plus  results,  we  believe 
that  the  errors  in  the  d-bands,  which  are  less  than  15!S,  are  a  result  of 
relaxation  effects,  which  would  be  expected  to  shift  the  d-bands  up 
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relative  to  the  top  of  the  valence  bands. 

In  Table  2c,  the  ionization  energies  for  the  cation  d-levels  of  GaP 
and  GaAs  are  compared  to  experimental  values.  This  information  is  not 
given  by  any  of  the  previous  band  calculations  mentioned .  The  differences 
between  theory  and  experiment  are  easily  accounted  for  by  relaxation  and 
surface  polarization  effects. 

In  Table  3.  the  lowest  transition  energies  are  compared  to  experiment. 
Clearly  the  Kartree-plus  conduction  bands  are  in  very  poor  agreement  v.’ith 
experiment.  To  test,  for  errors  in  the  treatment  of  the  self-energies,  v:e 
performed  a  calculation  for  A1P  with  the  self-energies  set  equal  to  zero. 
Tnis  calculation  yielded  conduction  bands  which  differ  from  those  with 

nonzero  self-energy  by  no  none  than  0.2  eV.  We  therefore  conclude  that  the 

errors  in  the  conduction  bands  arc  due  to  the  choice  of  the 

exchange-correlation  potential,  Vx„.  It  is  not  really  surprising  that  our 
choice  of  exchange-cornel atior.  potential  gives  good  results  for  the  valence 
bards  but  not  for  the  conduction  bards.  As  was  shown  by  Pantelides,  et 
£  1  •  ' ,  for  both  insulators  ar.d  wide  gap  reni conductors ,  correlation  effects 
in  valence  bands  are  a  result,  of  virtual  scattering  of  holes,  while- 
correlation  of  the  conduction  bands  comes  primarily  from  virtual  scattering 
of  electrons.  Since  the  mechanisms  for  correlation  in  valence  and 
conduction  bands  are  different,  it  is  unreasonable  to  expect  one 

exchange-correlation  potential  to  be  valid  for  both  cases.  V'e  suggest  that 
in  future  Hartree-plus  calculations,  it  should  be  possible  to  obtain  better 
results  by  using  two  separate  potentials  for  the  valence  and  conduction 
stat.es. 

6 .  CO.h’Oi  PSION 

Based  upon  the  results  given  by  Kunz,  et  at.1,  for  the  rare  gases  ar.d 


,  for  the  rare  gases 


the  results  for  III-V  semiconductors  reported  here,  v/e  conclude  that  the 
Hartree-plus  method  should  give  good  results  for  the  valence  bands  of  both 
insulators  and  semiconductors  when  used  with  the  RBKS  screened  Kohn-Sham 
exchange.  This  method  not  only  gives  good  results  for  the  valence  band 
widths,  it  also  gives  reasonable  values  for  the  core  d  levels  and  the 
ionization  potentials.  Based  upon  the  calculations  given,  we  suggest  that 
it  might  be  possible  to  find  a  single  exchange-correlation  potential  that 
could  be  used  to  obtain  good  conduction  bands  for  insulators  and 
semiconductors  in  a  two-potential  band  calculation.  Such  a  calculation 
would  use  one  potential  for  valence  states  and  a  second  potential  for 
conduction  states.  This  type  of  calculation  not  only  would  be  expected  to 
give  better  Hartree-plus  band  structures,  but  would  also  be  more  consistent 
with  what  we  know  about  correlation  in  semiconductors  and  insulators. 
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CHAPTER  1 


INTRODUCTION 


The  electronic  structure  of  surfaces  is  important  in  many 
scientific  and  technological  fields,  such  as  corrosion,  semiconductor 
devices,  and  catalysis.  However,  the  understanding  of  fundamental  solid 
surface  properties  has  lagged  behind  the  progress  in  understanding  of 
bulk  properties.  This  situation  has  been  a  result  of  the  lack  of 
experimental  methods  for  the  quantitative  characterization  of  the 
surface  structure,  and  the  theoretical  complexities  resulting  from  both 
the  loss  of  symmetry  at  the  surface  and  any  possible  reconstruction 
occuring  at  the  surface.  Recently  developed  experimental  techniques 
such  as  electron  energy  loss  spectroscopy  (ELS)  (  ),  two  photon 
picosecond  spectroscopy  (  ),  surface  extended  x-ray  absorption  fine 
structure  (EXAFS)  (  ),  and  extended  appearance  potential  fine 
structure  (EAPFS)  (  )  have  made  accessible  the  experimental  study  of 
the  structural  and  electronic  properties  of  real  surfaces.  The 
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development  of  powerful  computational  methods  has  made  possible  the 
theoretical  study  of  these  systems. 


1.1  THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE  ALUMINUM  (100)  SURFACE 

Aluminum  surfaces  oxidize  readily  upon  exposure  to  molecular 
oxygen.  Several  experimental  studies  have  been  performed  on  the  initial 
oxidation  of  aluminum  surfaces  (  ).  Ultraviolet  photoemission 
spectroscopy  (UPS)  (  )  and  Auger  spectroscopy  (  )  studies  have  been 
reported  for  many  crystal  orientations  of  aluminum,  including  the  (100) 
surface.  LEED  (  ),  EXAFS  (  ),  and  EAPFS  (  )  studies  have  been  reported. 

Ultraviolet  photoemission  experiments  (  )  have  reported  a  valence 
band  resonance  at  1.5  eV  below  the  Fermi  level  for  low  oxygen  exposure. 
This  resonance  is  attributed  to  the  oxygen  2p  band.  The  shift  of  the 
aluminum  2p  core  level  has  also  been  measured  for  both  low  and  high 
oxygen  exposures.  This  shift  toward  lower  energy  is  due  to  the  dipole 
moment  arising  from  the  charge  transfer  from  the  aluminum  substrate  to 
the  adsorbed  oxygen.  V/ork  function  measurements  have  also  been  reported 
(  ),  with  the  work  function  found  to  decrease  with  increasing  oxygen 
exposure.  This  decrease  has  been  interpreted  as  being  a  result  of 
incorporation  of  the  electronegative  oxygen  adatoms  beneath  the  aluminum 


surface. 
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LEED,  EXAFS,  and  LAPFS  techniques  are  able  to  probe  the  surface 
geo.T.etry  of  the  aluminum  plus  oxygen  adsorbate  system,  and  thus  these 
studies  are  of  special  interest  as  they  will  help  to  differentiate 
between  various  theoretical  models.  For  low  oxygen  coverage,  the 
reported  Al-0  bond  length  is  1.98  A  (  ),  in  agreement  with  previous 
theoretical  calculations  (  ).  For  higher  exposures,  a  value  of  1.88  A 
is  reported  (  ) ,  which  correspondsto  the  nearest-neighbor  distance  in 
AI2O3 . 

For  theoretical  studies  the  (100)  surface  is  ideal:  calculated 
electronic  properties  can  be  compared  to  existing  Auger  and  UPS  results, 
and  calculated  bond  positions  and  lengths  can  be  compared  with  LEED, 
EXAFS,  and  EAPF3  results. 

Lang  and  Williams  (4)  have  studied  theoretically  the  adsorption  of 
oxygen  onto  a  jellium  surface.  While  jellium  is  a  good  approximation  to 
bulk  aluminum,  it  is  not  adequate  to  describe  the  local  adsorption 
process.  For  this  reason,  cluster  studies,  which  are  particularly  well 
suited  to  describe  such  local  phenomena  (  ),  have  been  performed. 
Harris  and  Painter  (  )  and  Messmer  and  Salahub  (  ) ,  have  studied  the 
adsorption  of  atomic  oxygen  onto  the  hole  site  of  the  A1  (100)  surface, 
reporting  electronic  structures  which  agree  with  existing  UPS  data  for 
the  initial  adsorption  of  oxygen  onto  this  alimunum  surface.  Studies  of 
the  adsorption  of  oxygen  at  the  remaining  two  high-symmetry  sites  of 
this  surface  have  not  previously  been  reported.  An  accepted  explanation 
of  the  incorporation  of  oxygen  atoms  beneath  the  aluminum  surface  does 
not  yet  exist:  the  calculated  electronic  structure  and  binding 


potential  curves  for  the  on-top  and  bridge  sites  presented  here  should 
shed  light  on  this  problem. 

In  this  study,  the  interaction  between  adsorbate  atoms  (in  this 
case  oxygen)  and  the  aluminum  (100)  surface  is  theoretically 
investigated  by  modeling  the  surface  with  a  small  number  of  atoms  and 
then  using  the  unrestricted  Hartree-Fock  approximation  to  calculate  the 
electronic  structure  and  potential  energy  curves.  From  these  potential 
energy  curves  the  binding  energy,  bond  distance,  and  force  constants  of 
the  interaction  are  determined.  The  calculated  one-electron  orbitals 
give  added  insight  into  the  chemisorption  process,  and  calculated 
densities  of  states  are  compared  with  photoemission  data.  This 
calculational  procedure  is  able  to  accurately  describe  the  localized 
bonding  of  an  adsorbate  onto  a  surface  site,  the  aspect  of  chemisorption 
considered  most  important  (  ). 

The  localized  nature  of  the  chemisorptive  bond  is  indicated  by 
several  kinds  of  experimental  evidence.  An  adsorbed  atom  or  molecule 
can  be  observed  to  hop  from  one  localized  site  to  another  (  ). 
Infrared  spectroscopy  studies  have  found  that  the  vibrational  spectra  of 
intermediates  on  surfaces  are  often  very  similar  to  the  spectra  of 
isolated  molecules  (  ) .  These  experimental  results  reinforce  the 
intuitive  notion  that  the  chemisorptive  bond  is  similar  in  nature  to  the 
familiar,  localized  chemical  bond. 
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In  general,  the  loss  of  symmetry  in  the  direction  perpendicular  to 
the  surface  leads  to  the  choice  between  two  drastic  approximations  in 
theoretically  modeling  the  surface-adsorbate  interaction.  The 
maintenance  of  translational  symmetry  in  the  plane  parallel  to  the 
surface  greatly  simplifies  the  computational  difficulties  and  allows  for 
the  use  of  the  powerful  methods  developed  for  bulk  band  structure 
calculations.  The  alternative  approximation,  simulating  a  small  section 
of  the  surface  by  a  cluster  of  atoms,  can  more  effectively  describe  the 
localized  interaction  of  the  surface  with  a  single  atom. 

Numerous  calculations  of  the  electronic  structure  of  surfaces  using 
the  approximation  of  a  semi-infinite  surface  have  been  reported  (  ). 
The  major  drawback  of  this  approach  is  that,  in  general,  a  semi-infinite 
monolayer  of  adsorbed  atoms  must  be  considered  in  a  chemisorption 
calculation.  This  limitation  is  a  direct  result  of  the  symmetry 
assumed,  and  therefore  only  qualitative  calculations  can  treat  a  single 
adatom  on  a  semi-infinite  surface. 

By  relaxing  the  symmetry  requirements  of  the  semi-infinite  surface, 
one  gains  increased  flexibility  and  freedom.  For  instance,  atoms  can  be 
brought  down  over  a  variety  of  possible  bonding  sites.  Roughened  and 
stepped  surfaces,  exposed  corners,  and  small  particles,  which 
experimental  work  indicates  are  often  sites  favored  for  chemisorption, 
can  easily  be  simulated  by  the  cluster  approach.  Use  of  the 
semi-infinite  surface,  by  contrast,  is  limited  to  studies  of  the  perfect 


surface 


( 
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Once  the  appropriate  cluster  has  been  chosen  to  simulate  the 
surface,  the  electronic  structure  of  the  cluster-adsorbate  system  must 
be  calculated.  The  method  used  here  is  the  unrestricted  Hartree-Fock 
approximation,  from  which  electronic  wavef unctions,  electron  energy 
eigenfunctions,  and  the  total  electronic  energy  are  calculated.  This 
calculation  is  repeated  at  varying  adsorbate-surface  separations,  and 
the  potential  energy  curve  can  be  generated  from  the  differences  in 
total  electronic  energy. 

The  Hartree-Fock  approximation  as  the  method  of  calculation  is 
attractive  due  to  two  important  advantages.  The  method  is  a  convenient 
first  step  towards  an  aja  initio  solution  of  the  meny-particle 
Schrodir.ger  equation  for  a  fernion  system,  and  gives  a  useful 
zeroth-order  viavefunctior.  for  a  perturbation  calculation  of  the 
many-body  problem.  Second,  it  is  the  most  advanced  calculational  method 
in  which  elementary  physical  intuition  is  applicable.  In  more  complete 
calculations,  there  is  no  longer  a  one-to-one  correspondence  between 
particles  and  one-particle  wavefunctions,  and  the  independent  particle 
approximation  is  no  longer  applicable. 


1.2  THE  ELECTKONIC  PROPERTIES  OF  THE  SILICA  SURFACE 


The  electronic  structure  of  crystalline  silica  (Si02)  is  of 
technological  interest  in,  for  example,  the  manufacture  of  solid  state 
electronic  devices.  Consequently,  there  has  been  considerable 
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experimental  effort  in  this  area.  In  a  recent  review  of  experimental 
and  theoretical  results,  Griscora  (  )  has  compiled  much  of  the  existing 
data  on  the  bulk  properties  of  silica.  The  band  gap  has  been  found  to 
be  about  8.9  eV,  with  a  valence  band  width  of  approximately  11  eV. 
Ultra-violet  photo  emission  spectroscopy  (UPS)  studies  by  Ibach  and  Rowe 
(  )  have  detected  no  occupied  surface  states  within  the  band  gap. 
Williams  (  )  has  studied  silica  surfaces  grown  on  crystalline  Si  and 
found  traps  at  2  eV  below  the  bulk  silica  conduction  band  edge.  These 
traps  act  like  Coulomb  centers  with  a  positive  charge.  It  is  with  the 
silica  surface  and  its  defects  that  the  current  investigation  is 
concerned. 

Several  recent  experimental  studies  have  investigated  surface 
states  of  silica  (  ).  Schwidtal  (  )  has  found  a  radiation-induced 
feature  at  SI  eV  on  the  high-energy  side  of  the  Si  L23  VV  (V=valence) 
transition  in  the  Auger  electron  spectrum  (AES),  and  suggests  that  it  is 
due  to  an  Si  VD  (D=defect)  transition.  He  further  suggests  that 
this  defect  is  the  E^  center,  which  is  a  dangling  singly-occupied  sp^ 
orbital  of  silicon.  Ibach  and  Rowe  (  ),  Fujiwara  and  Ogata  (  ), 
Lieske  and  Hezel  (  ) ,  and  Bermudez  and  Ritz  (  )  have  measured  electron 
energy  loss  spectra  (ELS)  and  found  peaks  in  the  second  derivative 
spectrum  at  about  3.5,  5.0,  and  7.4  eV  (Fig.  1).  Ibach  and  Rowe  have 
attributed  these  transitions  to  a  partially  oxidized  surface  region, 
SiOx,  where  1<X'2.  Fujiwara  and  Ogata  (  )  have  shown  the  states 
associated  with  these  transitions  to  be  at  the  surface,  and  concluded 
they  were  due  to  metastable  SiOx.  Lieske  and  Hezel  (  )  have  associated 
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the  final  states  with  Si-Si  bonds  found  in  SiO^.  Most  recently, 
Bermudez  and  Ritz  (  )  have  studied  this  transition,  and  found  the  final 
state  to  be  due  to  "a  chemically  stable  species  formed  after  rupture  of 
the  silica  network,"  and  attributed  these  properties  to  a  surface  Si=0 
double  bond.  For  this  bond  to  exist,  the  silica  surface  must  undergo 
reconstruction.  Evidence  for  such  reconstruction  from  low  energy 
electron  diffraction  (LEED)  measurements  is  due  to  Janossy  and  Kenyhard 
(  ).  Also,  Hochstrasser  and  Antonini  (  )  have  observed  what  they 
believe  to  be  recombination  luminescence  due  to  the  rearranging  of  the 
surface  Si-0  bonds  immediately  after  cleaving.  They  measured  the 
lifetime  of  this  luminescence  to  be  less  than  10-6  seconcts. 

Theoretical  studies  of  the  silica  bulk  electronic  properties  have 
been  performed  (  ),  and  the  calculated  densities  of  states  are  in 
good  agreement  with  experiment,  except  for  the  AES  peak  at  91  eV  and  the 
ELS  peaks  at  3-3,  5.0,  and  6.8  eV.  both  of  these  discrepancies  can  be 
accounted  foy'  by  allowing  for  the  existence  of  islands  of  silicon  in  the 
bulk  Si02.  Bennett  and  Roth  (  )  have  calculated,  using  the  Huckel 
approximation,  the  electronic  properties  of  many  clusters  which 
approximate  silica,  and  in  particular  gave  attention  to  defects.  These 
calculations  were  performed  in  the  bulk;  however,  the  clusters  were 
sufficiently  small  that  they  can  be  interpreted  in  terras  of  surface 
defects,  or  in  terras  of  a  surface  SiOx  region.  For  oxygen  deficient 
clusters  with  a  relaxation  of  the  silicon  atoms  near  the  defect,  the 
calculated  energy  levels  showed  transitions  at  approximately  the  same 
energy  values.  These  are,  however,  probably  due  to  small  clusters  (2-3 


atoms)  of  silicon  in  the  SiOx.  Laughlin,  et  al.  (  ),  using  a  tight 
binding  method,  have  calculated  similar  defects  and  have  also  found 
states  which  appear  in  the  gap,  at  about  5  eV  above  the  valence  band 
edge.  They  have  also  studied  the  e'  center,  and  found  a  level  lying 
just  below  the  bulk  conduction  band.  Pollmann  and  Pantelides  (  )  have 
performed  a  calculation  using  a  Green's  function  formulation  on  the 
ideal  terminated  cubic  (^> -cristobalite)  surface.  Using  an  admittedly 
crude  model  for  the  surface,  they  found  no  states  in  the  band  gap. 

The  purpose  of  the  current  investigation  is  to  study  the  silica 
surface,  giving  particular  attention  to  the  three  ELS  peaks  observed  at 
3.5,  5,  and  7  eV.  The  cluster  model  is  again  employed,  for  the  same 
reasons  given  in  section  1.1.  Of  the  three  models  proposed  by 
experimentalists  to  account  for  these  transitions,  the  only  one  which 
has  previously  been  calculated  theoretically  is  theSiOx  model.  In 
chapter  are  presented  calculations  of  the  one-electron  energy  levels 
for  the  remaining  two  models,  the  E’  center,  and  the  Si=0  surface 
double  bond.  A  comparison  is  then  made  between  the  three  models  and 


experiment. 


METHOD  OF  CALCULATION 


An  exact  solution  to  the  wave  equation  for  a  large  cluster  of  atoms 
is  not  at  present  attainable.  In  fact,  such  a  solution,  while  nice  in  a 
mathematical  sense,  would  be  difficult  to  interpret  physically.  The  ab 
initio  Hartree-Fock  theory  provides  an  approximate  method  of  solution  to 
determine  the  electronic  structure  of  a  large  cluster  of  atoms  which  is 
both  easily  handled  by  contemporary  computers  and  well  suited  for 
straightforward  interpretation. 

2.1  HARTREE-FOCK  THEORY 

Ab  iniiia  Hartree-Fock  theory  is  used  to  calculate  the  electronic 
wavefunetions  and  energies  for  the  finite  clusters  studied  in  this  work. 
The  solids  studied  contain  atoms  of  low  atomic  number;  therefore 
relativistic  effects  are  small  and  will  be  ignored.  Using  the 
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nonrelativist ic  Schrodinger  theory,  the  many-body  wavefunction  T,  which 
is  a  function  of  electronic  and  nuclear  space-spin  coordinates,  is 
determined  by  the  equation 


where  H  is  the  many-body  Hamiltonian 
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Here  we  have  used  atomic  units,  where^h',  me,  and  e=1,  the  unit  of  length 
is  the  bohr  (1  bohr^-  0.529A),  and  the  energy  is  given  in  Hartree  (1  Hy 
1^27.2  eV).  Upper  case  letters  refer  to  nuclear  properties:  ~Rj  is  the 
position  of  the  I^*1  nucleus,  and  Mj  is  its  mass.  Lower  case  letters 
refer  to  electronic  properties:  “rj  is  the  coordinate  of  the  ith 
electron.  The  four-vector  x^  denotes  both  the  spatial  coordinates  and 
the  spin  of  the  ith  electron. 

The  Schrodinger  equation  is  simplified  by  employing  the 
Born-Oppenheimer  approximation  (  )  to  separate  the  nuclear  and 
electronic  coordinates.  The  total  wavefunction  is  assumed  to  be 


Page 


separable: 

ia,/)-  (2a) 

The  nuclear  wavef unction  '[(X),  is  now  a  solution  of  the  equation 


while  the  electronic  wavefunction,  is  a  solution  of  the  equation 

{f  <•?,*)  --Pet)  f  o?J)  (j 

A  solution  of  this  electronic  Schrodinger  equation  depends  on  the 
nuclear  coordinates  R  and  in  turn  defines  a  potential  energy  function  of 
the  coordinates: 


which  determines  the  notion  of  the  nuclei.  The  problem  is  now  reduced 
to  solving  the  electronic  Schrodinger  equation  for  a  particular  set  of 
nuclear  coordinates.  This  equation  cannot  in  general  be  solved  exactly; 
one  therefore  takes  refuge  in  the  independent  particle  model  for  a 
method  of  obtaining  an  approximate  solution. 

This  approach  is  mathematically  equivalent  to  a  separation  of 
variables  technique.  The  many-electron  wavefunction,  which  is  a 
functional  of  the  space-spin  coordinates  of  all  of  the  electrons,  is 
assumed  to  be  a  product  of  orbitals  which  are  functions  of  one-electron 
space-spin  coordinates.  However,  the  Pauli  principle  requires  that  the 
many-electron  wavefunction  be  antisymmetric  under  the  exchange  of  any 
two  electrons.  Therefore,  the  total  electronic  wavefunction  is 
approximated  by  an  antisyrametrized  product  of  one-electron  orbitals: 
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where  ^-is  the  antisymmetrizer.  Since  the  Schrodinger  equation  is  still 
not  in  general  solvable  analytically,  the  variational  principle  is 
introduced  to  obtain  approximate  solutions.  The  variational  principle 
states  that  the  normalized  expectation  value  of  the  Hamiltonian  using 
approximate  wavefunctions  is  a  rigorous  upper  bound  to  the  exact  energy. 
The  approximate  wavefunction  can  then  be  varied  to  minimize  the  energy 
and  thereby  be  the  best  approximation  to  the  exact  wavefunction. 
Applying  the  variational  principle  to  the  electronic  Schrodinger 
equation,  and  using  an  antisymmetrized  product  of  one-electron  orbitals, 
we  obtain  the  Hartree-Fock  equations: 


where 


l,-.  f;‘0 

1 


the  rock  operator: 


The  matrix  elements  are  Lagrange  multipliers,  and  have  been 
introduced  to  insure  the  orthonormality  of  the  one-electron  orbitals. 
In  general,  a  unitary  transformation  can  be  performed  on  the  orbital 
space  which  diagonalizes  this  matrix.  The  Hartree-Fock  equations  can 
then  be  written  in  the  standard  form: 


These  integro-differential  equations  are  normally  solved  via  an 
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iterative  procedure.  Thus,  the  electron  orbitals  are  determined  in  the 
average  field  of  all  the  other  electrons.  The  only  constraints  imposed 
on  the  wavefunction  are  that  the  orbitals  be  orthonormal  and  be 
functions  of  the  space-spin  coordinates  of  only  one  electron.  This 
level  of  approximation  is  called  generalized  Hartree-Fock  (GHF)  and  has 
yet  to  be  solved,  so  further  constraints  are  imposed  to  simplify  the 
calculation.  The  simplest  constraint  which  can  be  imposed  is  to  require 
that  the  orbitals  be  eigenfunctions  of  Sz,  i.e.,  the  orbitals  be  of  the 
form: 


ft  <?)  =  ; 
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where  < *  C r )  is  the  spatial  part  of  the  orbital,  and  u.  and  /?  are  the 
spin-up  and  spin-down  eigenfunctions  of  Sz,  respectively.  The  form  of 
the  Hartree-Fock  equations  is  unaffected  by  this  constraint,  and  this 
level  of  approximation  is  known  as  the  unrestricted  Hartree-Fock  (UHF). 
The  UHF  method  has  been  used  for  numerous  practical  calculations  of 
atoms,  molecules,  and  solids  (  ),  and  is  the  method  used  in  this 
study . 

Additional  restrictions  may  be  placed  upon  the  one-electron 
orbitals  in  order  to  further  simplify  the  calculations.  The  spatial 
parts  of  the  orbitals  may  be  required  to  be  symmetry  eigenfunctions  of 
the  cluster  under  consideration,  i.e.,  each  orbital  must  transform  as  an 
irreducible  representation  of  the  space  group  of  the  cluster.  Finally, 
for  doubly  occupied  orbitals  the  spin-up  and  spin-down  spatial  orbitals 
may  be  assumed  to  be  equivalent.  This  level  of  approximation  is  called 
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the  restricted  Hartree-Fock  method  (RHF).  One  advantage  of  RHF  is  that 
the  one-electron  orbitals  are  eigenfunctions  of  both  spatial  symmetry 
and  spin.  However,  for  singly  occupied  orbitals  complications  arise  in 
that  the  Fock  operator  differs  for  open  and  closed  shells  of  the  same 
symmetry,  thus  not  allowing  the  single  determinant  to  properly  describe 
the  wavefunction.  The  UHF  method,  while  not  necessarily  obtaining  a 
wavefunction  which  is  an  eigenfunction  of  spin  and  spatial  symmetry, 
does  give  a  lower  (i.e.,  better)  eigenfunction  of  energy.  Since  the 
"best”  solution  to  equation  2.6  is  the  wavefunction  which  gives  the 
lowest  energy,  the  UHF  method  has  therefore  been  used  in  this 
calculation. 


2.2  CLUSTER  MODEL 


Localized  phenomena  in  solids,  such  as  point  defects  or  surface 
adsorbates,  can  be  modeled  effectively  by  using  a  finite  cluster  of 
atoms  (  ).  Since  the  cluster  is  intended  to  represent  a  much  larger 
system,  one  must  apply  appropriate  boundary  conditions  for  the  cluster 
to  simulate  the  effects  of  the  environment.  These  have  been  discussed 
by  Kunz  and  Klein  (  ),  and  will  be  briefly  reviewed  here. 

The  solution  to  the  Hartree-Fock  equation  (2.11)  is  sought  for  a 
finite  cluster.  Suppose  the  region  of  the  cluster  is  termed  A,  and  the 
remainder  of  the  system,  the  environment  of  A,  is  termed  E.  The  problem 
is  to  partition  the  system  rigorously  into  a  cluster  and  an  environment. 
Let  us  consider  the  method  of  local  orbitals  of  Adams,  Gilbert,  and 


Page  fl y 

Kunz(  ).  Here  the  Fock  operator,  F,  of  the  entire  system  is  divided 
into  a  part  which  includes  the  kinetic  energy,  nuclear  attraction  of  the 
electrons  and  nuclei  inside  A,  along  with  the  electron-electron 
potential  for  the  electrons  assigned  to  A,  F^;  and  another  part  which 
is  the  potential  in  A  due  to  the  environment,  U^: 


It  is  desired  to  study  only  part  of  the  system,  so  instead  of  the  normal 
Hartree-Fock  equation  (2.11),  consider  instead  the  local  orbitals 
equation: 


where  W  is  an  arbitrary  Hermitian  operator.  Consider  now  an  ionic 
system.  The  potential  due  to  the  environment,  UA,  may  be  divided  into 

M 

two  parts:  V^1  is  an  ionic  (Madelung)  contribution  and  is  long  range, 

O  O 

and  is  the  remainder  and  is  short  range.  Let  W  =  V^.  Making  use 

of  the  projector  properties  of  ^  f, '  ,  and  solving  for  the  occupied 

orbitals,  equation  (2.14)  becomes 


Considering  only  the  orbitals  of  equation  (2.15)  which  lie  in  A,  the 
solutions  should  penetrate  only  weakly  into  E.  Because  does  not 
appreciably  penetrate  A,  and  because  ^  is  cancelled  by  j_  in  the 
limit  of  self-consistency,  the  approximate  equation  for  the  cluster, 
including  the  interaction  with  the  remainder  of  the  system,  is  then 
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This  approximate  equation  is  simply  the  equation  for  the  wavefunction  in 
the  cluster  of  interest  imbedded  in  the  Madelung  field  of  the  remainder 
of  the  solid.  The  eigenvalues'^  of  equation  (2.14)  or  the  approximate 
equation  (2.16)  represent  the  eigenvalues  C\  of  the  infinite  solid  only 
in  the  limit  that  the  orbitals  are  localized  in  the  region  A. 


A  basis  function  expansion  is  used  to  determine  the  unrestricted 
Hartree-Fock  orbitals.  Here  gaussian  type  orbitals  have  been  used  for 
this  expansion.  This  type  of  function  offers  the  advantage  of  having 
analytic  solutions  for  the  required  integrals;  however,  many  such 
functions  are  required  to  achieve  the  desired  accuracy.  Cartesian 
gaussian  functions  have  been  used.  The  spatial  part  of  these  functions 
has  the  form: 
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The  determination  of  the  gaussian  exponents,  for  atoms  and  ions  has 
been  the  subject  of  considerable  study,  and  tabulated  sets  exist  in  the 
literature.  An  all-electron  calculation  including  the  core  electrons  of 
aluminum  and  silicon  would  be  impractical,  so  the  ab  initio  effective 
potential  of  Topiol,  et  al.  (  )  is  used  to  replace  these  core 
electrons.  The  integrals  are  performed  using  the  standard  POLYATOM 
integrals  program,  and  the  unrestricted  Hartree-Fock  calculations  are 
done  with  the  G.  T.  Surratt  program  UHFONE. 
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THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE 
ALUMINUM  (100)  SURFACE 


Aluminum  crystalizes  in  the  FCC  form,  with  a  lattice  constant  of 
4.05  A.  Aluminum  is  a  metal,  and  the  bulk  electronic  structure  is 
approximated  very  well  by  the  "jellium"  model.  However,  attempts  to 
enploy  this  model  to  describe  the  local  chemisorption  process  (10)  have 
produced  results  which  are  not  consistent  with  experiment  (6).  Cluster 
techniques  have  been  used  to  study  the  chemisorption  process  (11,12,40), 
with  results  which  agree  well  with  experiment  if  care  is  taken  in 
choosing  the  cluster. 

Chemisorption  is  a  localized  process,  and  as  such  the  cluster  model 
should  provide  an  accurate  description.  The  difficulties  that  one 
encounters  here  are  resultant  from  the  fact  that  aluminum  is  a  metal, 
and  as  such  is  not  well  described  by  a  small  cluster  of  atoms  (fig. 
3.1).  One  must,  therefore,  take  into  account  the  size  of  the  cluster 
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used  to  describe  the  substrate.  These  effects  will  be  discussed,  along 
with  the  binding  energy  of  the  oxygen  atom  for  each  surface  site 
considered.  Calculated  one-electron  energy  levels  for  these  binding 
sites  are  compared  with  the  results  of  photoemission  experiments. 


3.1  RESULTS  OF  CALCULATIONS 


The  positions  of  the  aluminum  centers  in  the  clusters  are  chosen  to 
reproduce  the  unreconstructed  structure  of  the  aluminum  (100)  surface. 
It  has  been  assumed  here  that  the  lattice  does  not  undergo  signifigant 
distortions  at  the  surface.  This  has  been  shown  to  be  true  from 
comparison  of  experimental  and  theoretical  LEED  studies  of  clean 
aluminum  surfaces  (44). 

As  possible  surface  sites  for  the  adsorbed  oxygen  atom,  considered 
here  are  the  three  high-symmetry  points  of  the  (100)  surface  (fig. 
3.2).  These  are  a  position  directly  above  an  aluminum  center  of  the 
first  substrate  layer  (on-top  position);  the  position  central  to  four 
aluminum  atoms  of  the  first  layer,  above  a  second  layer  atom  (hole 
position);  and  above  (or  below)  the  midpoint  of  a  line  connecting  two 
nearest  neighbor  aluminum  centers  (bridge  position). 

The  interaction  of  oxygen  with  the  aluminum  surface  at  the  on-top 
position  has  been  modeled  here  using  two  clusters:  A10,  with  a  single 
aluminum  atom  representing  the  surface;  and  AlgO,  with  five  first  layer 
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atoms  and  four  second  layer  atoms  modeling  the  surface.  The  binding 
energy  of  the  oxygen  atom  to  the  aluminum  surface  is  shown  as  a  function 
of  its  distance  above  the  surface  in  figure  3-3*  For  both  surface 
clusters,  the  equilibrium  distance  is  about  3*5  bohr.  The  binding 
energies  differ  between  the  two  clusters;  this  energy  is  2.25  eV  for 
the  A10  cluster,  and  1.36  eV  for  the  AlgO  cluster.  The  charge  transfer 
from  the  aluminum  surface  to  the  oxygen  atom  is  about  0.7  electron  in 
both  cases. 

For  the  hole  site,  the  aluminum  surface  is  modeled  with  two 
clusters:  Al^o,  which  has  four  aluminum  centers  in  the  first  surface 
layer,  and  the  central  aluminum  in  the  second  layer;  and  AlgO,  which  is 
the  Al^O  cluster,  to  v.'hich  four  additional  second  layer  aluminum  centers 
are  added.  The  importance  of  including  the  second  layer  atoms  at  this 
site  has  been  shown  by  previous  theoretical  studies  of  the  interaction 
of  oxygen  with  metal  surfaces  (45).  Tne  oxygen  binding  energies  as  a 
function  of  the  distance  above  the  surface  plane  are  shown  for  these 
clusters  in  figure  3-4.  Both  clusters  have  a  minimum  total  energy 
(maximum  binding  energy)  when  the  oxygen  center  is  about  0.2  bohr  below 
the  plane  of  the  first  surface  layer.  The  cluster  with  only  one  second 
layer  aluminum  atom  is  found  to  bind  the  oxygen  adsorbate  by  5.4  eV, 
with  a  charge  transfer  to  the  oxygen  of  1.4  electrons.  The  larger 
cluster  is  found  to  bind  the  oxygen  by  4.8  eV,  with  a  charge  transfer  of 
1.3  electrons  to  the  adsorbed  oxygen. 
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The  binding  energies  computed  for  this  site,  as  well  as  those  for 
the  bridge  site,  must  be  carefully  extracted  from  the  calculated  total 
energies.  Since  the  oxygen  atom  penetrates  into  the  cluster  of  aluminum 
centers,  the  basis  functions  used  to  describe  the  electronic 
wavefunctions  associated  with  the  oxygen  adorbate  may  improve  the 
description  of  the  substrate  electronic  structure,  thus  lowering  its 
total  energy.  This  improvement  in  the  substrate  wavefunction  can  be 
integrated  into  the  caculation  of  the  binding  energy  by  simply  including 
the  basis  functions  of  the  oxygen  atom  when  calculating  the  total  energy 
of  the  substrate  cluster. 

The  bridge  site  has  been  modeled  here  with  three  different 
clusters.  The  first  cluster  consists  of  two  nearest  neighbor  aluminum 
centers  of  the  surface,  anc  an  oxygen  atom  above  the  midpoint  of  the 
line  connecting  these  centers.  The  second  cluster  contains  these  two 
centers,  as  well  as  two  additional  centers  beneath  them  from  the  second 
aluminum  layer.  For  the  final  cluster,  two  aluminum  centers  from  the 
third  layer  are  added.  The  inclusion  of  the  deeper  layer  aluminum 
centers  at  this  bonding  site  is  necessary  because  the  oxygen  atom 
penetrates  the  surface,  and  can  move  vertically  in  a  "tunnel"  between 
pairs  of  aluminum  centers  (see  fig.  3.2).  The  binding  energies  of  an 
oxygen  atom  interacting  with  these  aluminum  clusters  are  shown  in  figure 
3.5.  For  the  cluster  with  two  aluminum  centers,  the  oxygen  adsorbate  is 
bound  1.3  bohr  above  the  plane  of  the  surface.  A  charge  of  1.25 
electrons  is  transfered  to  the  oxygen  atom,  which  is  bound  by  3-1  eV. 
With  the  inclusion  of  the  two  second  layer  aluminum  centers  in  the 
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cluster,  the  equilibrium  position  of  the  oxygen  center  moves  to  1.65 
bohr  beneath  the  surface,  which  is  just  above  halfway  between  the  two 
aluminum  layers.  Tne  binding  energy  is  found  to  be  3.8  eV,  with  a 
charge  transfer  of  I.1!  electrons  to  the  oxygen.  Inclusion  of  the  third 
aluminum  layer  produces  no  signifigant  change  in  these  results; 
however,  this  cluster  does  show  that  the  oxygen  center  does  not  want  to 
move  deeper  into  the  bulk  material,  as  the  total  energy  of  the  system  is 
found  to  be  lowest  for  an  oxygen  center  position  just  below  the  first 
surface  layer. 


3.2  DISCUSSION  OF  RESULTS 


For  the  initial  stages  of  adsorption,  the  oxygen  adsorbate  is 
clearly  bound  below  the  aluminum  surface.  The  binding  is  strongest  at 
the  hole  site.  An  electric  charge  of  about  1.35  electrons  is  also 
transfered  to  the  adsorbed  oxygen.  These  results  are  in  agreement  with 
previously  reported  work  function  studies,  in  which  a  reported  decrease 
in  the  work  function  of  the  aluminum  surface  with  increasing  oxygen 
exposure  (up  to  one  monolayer  of  coverage)  has  been  attributed  to  the 
penetration  of  negative  oxygen  ions  beneath  the  aluminum  surface  (6). 

Tne  ultra-violet  photoemission  spectrum  has  been  reported  for  low 
oxygen  coverage  on  the  aluminum  (100)  surface  (5);  this  spectrum  is 
shown  in  figure  3.6.  Also  shown  is  the  projected  density  of  states  for 
the  oxygen  adsorbate  previously  calculated  by  Messmer  and  Salahub  (12) 


for  an  oxygen  position  at  the  aluminum  surface  at  the  hole  site.  These 
authors  did  not  calculate  the  electronic  structure  for  an  oxygen 
position  below  the  surface,  nor  did  they  study  the  interaction  of  an 
oxygen  adsorbate  with  any  of  the  other  high-symmetry  (100)  surface 
sites. 

As  reported  elsewhere  (40),  the  one-electron  energy  levels  of  this 
calculation  are  in  agreement  with  experiment  and  previous  theory.  These 
energy  levels  are  shown  in  figure  3.6  for  both  the  hole  and  bridge 
sites,  with  the  oxygen  position  corresponding  to  the  minimum  total 
energy  for  the  cluster.  For  both  sites,  the  energy  levels  reported  here 
are  those  of  the  largest  cluster  calculated. 

Photoemission  experiments  studying  the  aluminum  core  levels  at  low 
oxygen  exposures  have  been  reported,  with  a  reported  shift  in  the 
aluminum  2p  energy  of  1.3eV  toward  lower  energy  (5).  In  the  present 
calculation,  the  shift  of  the  aluminm  core  levels  is  found  to  be  1.1  eV 
ir.  the  same  direction.  The  cause  of  this  core  level  shift  is  the 
electric  dipole  produced  by  the  charge  transfer  from  the  aluminum 
substrate  to  the  oxygen  atom  upon  adsorption.  It  has  been  noted  by 
Flodstrom,  et^  al.  that  the  aluminum  core  level  shift  increases  to  a 
value  of  2.6  eV  at  higher  oxygen  exposure  (about  one  monolayer  of 
coverage),  thus  lowering  these  orbital  energies  to  the  values  found  in 
AI2O3.  These  results  are  evidence  for  a  two  stage  oxidation  process, 
with  the  first  step  being  chemisorption  at  the  hole  site,  and  the  second 
step  being  incorporation  of  oxygen  beneath  the  aluminum  surface. 
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Hie  agreement  between  the  calculated  one-electron  energies  for  the 
hole  site  presented  here,  along  v.’ith  those  of  Messmer  and  Salahub,  and 
reported  experimental  UPS  spectra  tend  to  support  the  idea  that  the 
initial  chemisorption  of  oxygen  occurs  at  the  hole  site,  as  one  would 
expect  from  the  binding  energies  reported  in  section  3.1. 

The  hole  site  is  also  favored  by  extended  appearance  potential  fine 
structure  (EAPFS)  measurements.  Using  this  technique,  den  Boer,  et  al. 
(*J)  have  reported  a  nearest  neighbor  oxygen-aluminum  distance  of  1.98  A, 
or  3.74  bohr.  Messmer  and  Salahub  have  reported  a  value  of  2.02  A  (3.82 
bohr),  and  the  results  of  the  calculations  presented  here  place  the 
oxygen  at  1.92  A  (3.63  bohr)  from  the  nearest  aluminum  center  for  the 
hole  site.  For  higher  oxygen  exposure,  den  Boer^et  al.  have  reported 
an  experimental  value  of  this  nearest  neighbor  distance  of  1.91  A.  This 
figure  is  probably  an  average  of  the  two  nearest  neighbor  spacings  found 
in  A.I2O3 ;  these  being  1.86  A  (3.51  bohr)  and  1.97  A  (3*72  bohr).  The 
longer  of  these  spacings  corresponds  to  oxygen  binding  at  the  hole  site, 
while  the  shorter  corresponds  to  oxygen  adsorption  at  the  slightly  less 
energetically  favorable  bridge  site.  The  present  calculation  places  the 
nearest  neighbor  distance  at  the  bridge  site  at  1.71  A,  or  3.25  bohr. 

The  results  presented  above,  along  with  previous  theoretical  and 
experimental  results,  show  that  for  the  initial  interaction  of  oxygen 
with  the  aluminum  (100)  surface,  chemisorption  takes  place  at  the  hole 
site.  The  present  studies  have  extended  previous  theoretical 
understanding  of  the  interaction  of  oxygen  with  this  surface  to  include 
the  two  additional  high-symetry  (100)  sites,  l/ork  function,  UPS,  and 
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EAPFS  studies  have  reported  that  the  nature  of  this  interaction  changes 
when  the  oxygen  coverage  approaches  one  nonolayer;  the  theoretical 
results  presented  here  lead  one  to  conclude  that  at  this  coverage, 
oxygen  begins  to  adsorb  at  the  bridge  site,  since  the  energitically  more 
favorable  hole  sites  are  already  occupied.  The  chemisorption  of  oxygen 
at  the  bridge  site  allows  for  the  incorporation  of  the  oxygen  adatoms 
beneath  the  aluminum  surface  at  high  levels  of  oxygen  exposure,  leading 
to  the  formation  of  the  oxide  AI2O3. 


THE  ELECRONIC  PROPERTIES  OF  THE  SILICA  SURFACE 


Silica  occurs  in  six  different  crystalline  forms,  as  well  as  the 
amorphous  form.  All  forms  of  Si02  are  based  on  the  SiOjj  tetrahedron, 
and  differ  in  the  Si-O-Si  bonding  angle.  The  crystalline  form  most 
often  studied  experimentally  is  X-quartz,  and  it  is  this  form  which  we 
have  used  here,  t^-quartz  is  hexagonal  in  structure,  with  three  Si02 

m 

molecules  in  each  unit  cell.  The  Si-O-Si  bond  angle  is  IM  . 

The  silica  surface  may  be  described  accurately  using  a  cluster  of  a 
small  number  of  atoms.  Indeed,  calculations  on  the  bulk  material  (  ) 
have  used  the  cluster  approach.  To  take  into  account  the  effect  of  the 
neighboring  bulk,  appropriate  boundary  conditions  are  imposed,  as  has 
been  discussed  in  chapter  2.  In  the  case  of  silica,  the  cluster  is 
embedded  in  a  point  charge  array,  with  the  charge  transfer  determined  in 
a  self-consistent  manner.  This  gives  the  proper  electrostatic  potential 
in  the  region  of  the  surface  cluster,  and  also  provides  for  charge 
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neutrality.  Previous  calculations  have  either  employed  a  free  space 
termination,  or  imposed  periodic  boundary  conditions  which  require 
defects  to  be  present  in  all  unit  cells  of  the  crystal. 

The  clusters  used  here  for  the  bulk  calculations  were  an  SiO^ 
tetrahedron,  and  an  Si-O-Si  molecule.  Both  clusters  were  embedded  in  a 
point  charge  array  of  3*3*3  unit  cells  in  size.  The  calculations  of  the 
surface  were  done  using  a  cluster  consisting  of  the  surface  silicon 
atom,  its  two  neighboring  oxygen  atoms,  and  the  two  adjacent  silicon 
atoms  to  form  two  joined  Si-O-Si  molecules.  This  cluster  was  embedded 
in  a  point  charge  array  of  3*3*2  unit  cells.  The  position  of  the 
surface  silicon  atom  was  adjusted  in  order  to  minimize  the  total  energy. 
An  oxygen  atom  was  then  placed  above  the  surface  silicon  in  order  to 
study  the  surface  bonding  state. 

The  one-electron  energy  levels  calculated  for  the  bulk  SiC>2  are  in 
good  agreement  with  experiment  and  previous  theory  (fig.  4.1).  The 
calculated  band  gap  is  9.4  eV,  as  compared  to  8.9  eV  for  experiment. 
The  valence  band  width  is  calculated  to  be  9.2  eV,  while  experimental 
measurements  give  about  11  eV.  The  oxygen  2s  band  is  found  at  30.3  eV 
below  the  conduction  band,  while  experiment  places  it  28  eV  below.  The 
valence  bands  are  found  to  be  mostly  oxygen  2p,  with  some  silicon  3s  and 
3p  character  mixed  in.  This  is  expected,  and  agrees  with  the  charge 
transfer  of  nearly  one-half  of  an  electron  from  each  silicon  to  each  of 
4  oxygen  atoms,  or  a  configuration  of  Si^+o",  as  determined  by  a 
Kulliken  population  analysis.  This  agrees  with  the  observation  that 
silica  has  both  ionic  and  covalent  properties. 
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For  the  unreconstructed  surface,  the  E'  center  is  found  to  have  a 

s 

singly  occupied  energy  level  1.7  eV  above  the  top  of  the  bulk  valence 
band  (fig.  4.2).  Stephenson  and  Binkow3ki  (  )  have  observed,  usir.g 

XPS,  an  occupied  energy  level  at  0.75  eV  above  the  bulk  valence  band 
edge  for  samples  cleaved  in  vacuum.  They  believed  this  level  to  be  an 
intrinsic  feature  of  the  bulk  SiC^.  However,  their  sampling  depth  is 
only  about  30  A,  so  it  is  probable  that  this  level  is  an  E'  center  at  or 
near  the  surface  (  ).  The  method  of  sample  preparation  which  they  have 
used  (grinding)  could  easily  have  caused  these  broken  bonds  to  be  formed 
(  ).  It  is  energetically  favorable  for  an  oxygen  atom  to  bind  to  this 

surface  silicon  atom,  with  a  binding  energy  of  5.1  eV,  which  is 
considerably  greater  than  the  2.6  eV  per  oxygen  atom  necessary  to 
dissociate  molecular  oxygen.  We  shall  now  turn  our  attention  to  this 
configuration,  which  is  the  Si=0  surface  double  bond. 

The  occupied  one-electron  energy  levels  of  the  surface  state  are 

all  below  the  top  of  the  bulk  valence  band.  This  is  consistent  with 

optical  studies,  which  have  found  no  occupied  surface  states  in  the  band 

gap  (  ).  The  lowering  of  the  valence  band  can  be  attributed  to  the 

O-Si-O  bond  angle  at  the  surface  being  greater  than  the  perfect 
o 

tetrahedral  109  angle.  The  bottom  of  the  conduction  band  for  the 
surface  is  found  to  be  4.6  eV  below  the  bulk  conduction  band  edge.  The 
occupied  valence  surface  states  can  be  described  as  bonding  and 
non-bonding  states  between  the  adsorbed  oxygen  and  the  surface  silicon 
atom  (fig.  4.3a»b).  When  a  surface  electron  is  excited  out  of  the 
valence  band,  an  ariti-bonding  orbital  (fig.  4.3c)  drops  out  of  the 
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conduction  band  and  is  singly  occupied.  Its  energy  lowers  considerably 
due  to  relaxation,  and  is  found  at  I.1)  eV  above  the  top  of  the  bulk 
valence  band  (fig.  4.2).  The  position  of  the  oxygen  atom  is  0.1  a.  u. 
farther  out  from  the  surface  silicon  atom  for  this  excited  state, 
relative  to  its  ground  state  position.  This  energy  difference  is  in 
good  agreement  with  Bermudez  and  Ritz  (  ) ,  and  the  transition  from  a 
non-bonding  to  an  anti-bonding  orbital  of  the  surface  bond  is  as  they 
have  described.  The  energy  change  associated  with  this  transition  is 
6.1  eV,  which  is  in  agreement  with  the  experimental  value  of  5.0  eV  if 
correlation  effects  are  considered.  The  7.4  eV  transition  can  be 
understood  in  terms  of  valence  band  structure,  with  transitions  from 
levels  in  the  bonding  part  of  the  valence  band  to  the  same  anti-bonding 
orbital. 

It  is  noted  that  this  model,  the  Si=Q  surface  bond,  does  not 
predict  the  surface  electronic  transition  at  3-5  eV  that  is  seen  in  the 
second  derivative  ELS  spectra.  It  has  been  shown  by  Gallon  and 
Underwood  (  )  that  this  peak  is  in  fact  an  artifact  of  the  second 
derivative  mode  of  detection,  and  it  is  not  seen  in  the  non-derivative 
ELS  spectrum.  They  argue  that  the  3*5  eV  peak  is  produced  by  the 
overlap  of  the  "wings"  of  the  primary  and  the  5.0  eV  peaks.  The 
calculations  presented  here  agree  with  this  interpretation. 

One  remaining  point  of  controversy  is  the  difference  between  XPS 
and  UPS  in  discribing  the  top  of  the  valence  bands  of  SiOg  (fig.  4.4) 
There  exists  at  this  time  no  agreement  as  to  the  cause  of  this 
discrepancy.  The  calculations  presented  here  do,  however,  favor  one 
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explaination.  The  two  spectra  are  spatially  sensitive  to  different 
parts  of  the  sample:  XPS  is  sensitive  to  the  bulk  for  initial  states 
with  small  binding  energy,  while  40.8  eV  UPS  is  surface  sensitive.  In 
both  the  XPS  and  UPS  techniques,  the  kinetic  energy  of  the  photoelectron 
is  sufficiently  large  that  the  final  state  may  be  approximated  by  a 
plane  wave.  Therefore  the  difference  between  the  XPS  and  the  UPS 
spectra  must  reflect  the  difference  between  the  bulk  and  surface  valence 
band  densities  of  states.  Comparison  between  the  XPS  spectrum  and  the 
calculated  bulk  density  of  states,  and  between  the  UPS  spectrum  and  the 
calculated  surface  density  of  states  (for  the  reconstructed  surface) 
shows  that  such  an  explanation  does  indeed  account  for  the  difference  in 
the  two  spectra  (fig.  4.4). 

The  reconstructed  silica  surface  has  been  studied  in  order  to 
determine  the  nature  of  the  low  energy  loss  ELS  peaks.  We  have  found 
agreement  between  the  theoretically  calculated  one-electron  energy 
levels  presented  here  and  previous  experimental  results.  Of  the  three 
models  which  can  account  for  the  peaks,  namely  a  partly  oxidized  surface 
region  SiOx,  the  center,  and  the  Si=0  double  bond,  we  conclude  that 
the  last  is  most  likely.  Support  for  this  conclusion  comes  from 
experimental  evidence  for  reconstruction,  and  from  the  large  binding 
energy  of  the  oxygen  atom  to  the  surface. 
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CONCLUSIONS 


In  this  study,  the  sb  in litio  unrestricted  hartree-Fock  method  has 
been  used  to  calculate  the  bulk  and  surface  electronic  structure  of  of 
silica,  and  to  study  the  initial  interaction  of  oxygen  with  the  aluminum 
(100)  surface.  In  both  of  these  materials,  the  cluster  model  has  been 
shown  to  be  a  useful  tool  in  the  calculation  of  the  electronic 
structures. 

In  studying  the  interaction  af  oxygen  with  the  aluminum  (100) 
surface,  the  size  of  the  aluminum  cluster  used  had  a  considerable  effect 
on  the  results.  Since  aluminum  is  a  metal,  a  fairly  large  cluster  of 
atoms  is  needed  to  describe  the  substrate  accurately.  That  the  cluster 
model  worked  at  all  for  this  calculation  is  due  to  the  fact  that 


chemisorption  is  a  local  process 


( 
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It  has  been  shown  that  an  oxygen  adsorbate  is  bound  most  strongly 
at  the  hole  site  of  the  clean  aluminum  surface.  The  equilibrium 
position  of  the  adsorbed  atom  is  about  0.2  bohr  below  the  plane  of  the 
surface.  A  charge  of  1.3  electrons  is  trar.sfered  from  the  aluminum 
substrate  to  the  oxygen  adsorbate.  Thi'S  charge  transfer  creates  an 
electric  dipole  which  causes  the  aluminum  core  levels  to  be  shifted 
downward  in  energy  by  1.1  eV.  Previous  experiments  studying  the  work 
function,  UPS  spectrum,  and  the  EAPFS  have  reported  results  for  the 
initial  adsorption  of  oxygen  onto  this  surface  which  are  in  agreement 
with  the  results  of  the  calculations  presented  here. 

These  experimental  results  have  also  shown  that  the  nature  of  the 
adsorption  process  changes  at  about  one  monolayer  of  coverage.  The 
theoretical  studies  reported  here  lead  one  to  conclude  that  at  this 
coverage  the  energetically  most  favored  sites  for  chemisorption,  the 
hole  sites,  become  fully  occupied,  and  adsorption  continues  at  the 
bridge  site.  It  is  adsorption  of  oxygen  at  this  site  which  allows  for 
incorporation  of  oxygen  beneath  the  aluminum  surface,  and  leads  to  the 
formation  of  the  oxide  AI2O3. 

Also  reported  here  is  a  study  of  the  bulk  and  surface  electronic 
structure  of  silica  (e^-quartz).  The  cluster  model  was  again  used.  The 
bulk  solid  and  the  surface  were  both  modeled  with  clusters  of  a  few 
atoms,  plus  the  appropriate  boundary  conditions.  In  this  case,  the 
boundary  condition  imposed  was  a  point  ion  array  to  provide  for  charge 
neutrality  and  to  provide  the  correct  Madelung  field. 


m 

For  the  silica  bulk,  the  calculated  one-electron  energy  levels  were 
compared  with  experimental  XPS  spectra,  and  the  agreement  was  found  to 
be  excellent,  as  reported  in  chapter  >i.  For  the  reconstructed  silica 
surface,  the  structure  of  the  Si=0  double  bond  was  found  to  be  the 
energetically  most  favorable  surface  configuration.  The  electronic 
structure  of  this  surface  bond  is  able  to  explain  the  low  energy  ELS 
peaks  seen  experimentally .  This  calculations  reports  a  value  of  6.1  eV 
for  the  lowest  energy  transition,  as  compared  to  5.0  eV  for  the 
experimental  result. 

Finally,  the  difference  in  the  reported  valence  band  structures  as 
measured  in  XPS  and  UPS  experiments  is  explained.  Tne  XPS  method  is 
sensitive  to  the  bulk  material,  while  UPS  is  surface  sensitive. 
Comparison  of  these  spectra  with  the  calculated  valence  band  energy 
levels  shows  good  agreement  between  the  theory  presented  here  and  the 
experimenal  spectra. 

The  ab  initio  unrestricted  Hartree-Fock  method,  along  with  the 
cluster  model,  has  been  shown  to  accurately  describe  the  local 
electronic  properties  cf  many  systems.  In  this  report,  this  method  has 
been  shown  to  be  successful  in  describing  the  adsorption  of  oxygen  onto 
the  aluminum  surface,  and  in  describing  the  bulk  and  surface  electronic 
structure  of  SiC>2.  This  technique  has  also  been  used  to  successfully 
calculate  the  electronic  structure  of  semiconducting  polymers,  including 
defects  and  impurities  (%),  and  in  the  description  of  the  excitonic 
structure-  of  crystalline  silicon  (*I7). 
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